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Gravity solutions dual to d- dimensional field theories at finite charge density have a near-horizon 
region which is AdS2 x M — 1 . The scale invariance of the AdS? region implies that at low energies 
the dual field theory exhibits emergent quantum critical behavior controlled by a (0+ 1) -dimensional 
CFT. This interpretation sheds light on recently-discovered holographic descriptions of Fermi sur- 
faces, allowing an analytic understanding of their low-energy excitations. For example, the scaling 
behavior near the Fermi surfaces is determined by conformal dimensions in the emergent IR CFT. 
In particular, when the operator is marginal in the IR CFT, the corresponding spectral function is 
precisely of the "Marginal Fermi Liquid" form, postulated to describe the optimally doped cuprates. 
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I. INTRODUCTION 

The AdS/CFT correspondence [1] has opened new av- 
enues for studying strongly-coupled many-body phenom- 
ena by relating certain interacting quantum field theories 
to classical gravity (or string) systems. Compared to con- 
ventional methods of dealing with many-body problems, 
this approach has some remarkable features which make 
it particularly valuable: 

1. Putting the boundary theory at finite temperature 
and finite density corresponds to putting a black 
hole in the bulk geometry. Questions about com- 
plicated many-body phenomena at strong coupling 
are now mapped to single- or few-body classical 
problems in a black hole background. 
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2. Highly dynamical, strong-coupling phenomena in 
the dual field theories can often be understood on 
the gravity side using simple geometric pictures. 

3. At small curvature and low energies a gravity the- 
ory reduces to a universal sector: classical Einstein 
gravity plus matter fields. Through the duality, this 
limit typically translates into the strong-coupling 
and large-iV limit of the boundary theory, where 
iV characterizes the number of species. Thus by 
working with Einstein gravity (plus various matter 
fields) one can extract certain universal properties 
of a large number of strongly coupled quantum field 
theories. 

In this paper following [2, 3] we continue the study of non- 
Fermi liquids using the AdS/CFT correspondence (see 
also [4]). 

Consider a (i-dimensional conformal field theory (CFT) 
with a global U(l) symmetry that has an AdS gravity 
dual. Examples of such theories include the Af = 4 super- 
Yang-Mills (SYM) theory in d = 4, the Af = 8 M2 brane 
theory in d = 3, the (2,0) multiple M5-brane theory in 
d = 6, and many others with less supersymmetry. With 
the help of the AdS/CFT correspondence, many impor- 
tant insights have been obtained into strongly coupled 
dynamics in these systems, both near the vacuum and 
at a finite temperature. In particular, as a relative of 
QCD, thermal Af = 4 SYM theory has been used as 
a valuable guide for understanding the strongly coupled 
Quark-Gluon Plasma of QCD. 

It is also natural to ask what happens to the result- 
ing many-body system when we put such a theory at a 
finite U(l) charge density (and zero temperature). Im- 
mediate questions include: What kind of quantum liquid 
is it? Does the system have a Fermi surface? If yes, is it 
a Landau Fermi liquid? A precise understanding of the 
ground states of these finite density systems at strong 
coupling should help expand the horizon of our knowl- 
edge of quantum liquids, and may find applications to 
real condensed matter systems. 

On the gravity side, such a finite density system is 
described by an extremal charged black hole in d + 1- 
dimensional anti-dc Sitter spacetime (AdSd+i) [5]. The 
metric of the extremal black hole has two interesting 
features which give some clues regarding the nature of 
the system. The first is that the black hole has a fi- 
nite horizon area at zero temperature, suggesting a large 
ground state degeneracy (or approximate degeneracy) in 
the large ./V limit. The second is that the near hori- 
zon geometry is given by AdS2 x R d_1 , which appears 
to indicate that at low frequencies the boundary system 
should develop an enhanced symmetry group including 
scaling invariance. In particular, it is natural to expect 
that quantum gravity (or string theory) in this region 
may be described by a boundary CFT. It has been ar- 
gued in [6] that the asymptotic symmetry group of the 
near horizon AdS2 region is generated by a single copy 
of Virasoro algebra with a nontrivial central charge, sug- 



gesting a possible description in terms of some chiral 2d 
CFT . 

More clues to the system were found in [3] from study- 
ing spectral functions of a family of spinor operators (fol- 
lowing earlier work of [2]): 

1. The system possesses sharp quasi-particle-likc 
fcrmionic excitations at low energies near some dis- 
crete shells in momentum space, which strongly 
suggests the presence of Fermi surfaces. 1 In par- 
ticular, the excitations exhibit scaling behavior as 
a Fermi surface is approached with scaling expo- 
nents different from that of a Landau Fermi- liquid 2 . 
The scaling behavior is consistent with the general 
behavior discussed by Scnthil in [7] for a critical 
Fermi surface at the critical point of a continuous 
metal-insulator transition. 

2. For a finite range of momenta, the spectral func- 
tion becomes periodic in log uj in the low frequency 
limit. Such log-periodic behavior gives rise to a dis- 
crete scaling symmetry which is typical of a com- 
plex scaling exponent. 

Note that the above scaling behavior is emergent, a con- 
sequence of collective dynamics of many particles, not re- 
lated to the conformal invariance of the UV theory which 
is broken by finite density. 

The results of [3] were obtained by solving numerically 
the Dirac equation for bulk spinor fields dual to boundary 
operators and it was not possible to identify the specific 
geometric feature of the black hole which is responsible 
for the emergence of the scaling behavior. Nevertheless, 
as speculated in [3] , it is natural to suspect that the AdS2 
region of the black hole may be responsible. 

In this paper, we show that at low frequencies 3 , re- 
tarded Green functions 4 of generic operators in the 
boundary theory exhibit quantum critical behavior. This 
critical behavior is determined by the AdS2 region of the 
black hole; assuming it exists, a CFTi dual to this region 
of the geometry (which we will call the 'IR CFT') can be 
said to govern the critical behavior. 



1 Note that since the underlying system is spherically symmetric, 
the Fermi surfaces are round. 

2 In [4] a different family of operators were studied at finite tem- 
perature. The authors concluded there that the scaling behavior 
resembles that of a Fermi liquid. We will discuss those operators 
in section VI D. 

3 By "low frequency" we mean frequency close to the chemical 
potential. How this arises from the AdS/CFT dictionary will be 
explained below. 

4 We focus on retarded Green function as its imaginary part di- 
rectly gives the spectral function which reflects the density of 
states which couple to an operator. It is also the simplest ob- 
servable to compute in the Lorentzian signature in AdS/CFT. 
The scaling behavior is of course also present in other types of 
correlation functions. We also expect similar scaling behavior to 
exist in higher-point functions which will be left for future study. 
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The spirit of the discussion of this paper will be similar 
to that of [3]; we will not restrict to any specific theory. 
Since Einstein gravity coupled to matter fields captures 
universal features of a large class of field theories with 
a gravity dual, we will simply work with this universal 
sector, essentially scanning many possible CFTs. 5 

The role played by the IR CFT in determining the 
low-frequency form of the Green's functions of the d- 
dimcnsional theory requires some explanation. Each op- 
erator O in the UV theory gives rise to a tower of op- 
erators 0£ in the IR CFT labeled by spatial momentum 

k. The small 10 expansion of the retarded Green function 
Gu(uj,k) for O contains an analytic part which is gov- 
erned by the UV physics and a non-analytic part which 
is proportional to the retarded Green function of 0£ in 
the IR CFT. What kind of low-energy behavior occurs 
depends on the dimension 5k of the operator in the 
IR CFT and the behavior of Gr(lu = 0, k). e For example, 
when 8k is complex one finds the log-periodic behavior 
described earlier. When Gr(uj — 0, k) has a pole at some 
finite momentum \k\ = kp (with Sk F real), one then finds 
gapless excitations around \k\ = kp indicative of a Fermi 
surface. 

Our discussion is general and should be applicable to 
operators of any spin. In particular both types of scal- 
ing behavior mentioned earlier for spinors also applies 
to scalars. But due to Bose statistics of the operator in 
the boundary theory, this behavior is associated with in- 
stabilities of the ground state. In contrast, there is no 
instability for spinors even when the dimension is com- 
plex. 

Our results give a nice understanding of the low-energy 
scaling behavior around the Fermi surface. The scaling 
exponents are controlled by the dimension of the corre- 
sponding operator in the IR CFT. When the operator is 
relevant (in the IR CFT), the quasi-particle is unstable. 
Its width is linearly proportional to its energy and the 
quasi-particle residue vanishes approaching the fermi sur- 
face. When the operator is irrelevant, the quasi-particle 
becomes stable, scaling toward the Fermi surface with 
a nonzero quasi-particle residue. When the operator is 
marginal the spectral function then has the form for a 
"marginal Fermi liquid" introduced in the phenomcno- 
logical study of the normal state of high T c cuprates [8] . 

It is also worth emphasizing two important features of 
our system. The first is that in the IR, the theory has not 
only an emergent scaling symmetry but an SL(2, R) con- 
formal symmetry (maybe even Virasoro algebra). The 
other is the critical behavior (including around the Fermi 
surfaces) only appears in the frequency, not in the spatial 



momentum directions. 

The plan of the paper is as follows. In §2, we introduce 
the charged AdS black hole and its AdS2 near-horizon 
region. In §3, we determine the low energy behavior of 
Green's functions in the dual field theory, using scalars as 
illustration. The discussion for spinors is rather parallel 
and presented in Appendix A. In §IV-VI, we apply this 
result to demonstrate three forms of emergent quantum 
critical behavior in the dual field theory: scaling behavior 
of the spectral density (§IV) , periodic behavior in log ui at 
small momentum (§V), and finally (§VT) the Fermi sur- 
faces found in [3] . We conclude in § VII with a discussion 
of various results and possible future generalizations. We 
have included various technical appendices. In particu- 
lar in Appendix D we give retarded functions of charged 
scalars and spinors in the AdS2/CFTi correspondence. 



II. CHARGED BLACK HOLES IN ADS AND 
EMERGENT INFRARED CFT 

A. Black hole geometry 

Consider a ci-dimcnsional conformal field theory (CFT) 
with a global U{1) symmetry that has a gravity dual. At 
finite charge density, the system can be described by a 
charged black hole in d + 1-dimcnsional anti-de Sitter 
spacetime (AdS^+i) [5] with the current J M in the CFT 
mapped to a U(l) gauge field Am in AdS. 

The action for a vector field Am coupled to AdS^+i 
gravity can be written as 



S 



d(d-l) R 



R 2 



q-FmnF 

9 2 f 



AIN 



(1) 

where gp is an effective dimensionless gauge coupling 7 
and R is the curvature radius of AdS. The equations of 
motion following from (1) are solved by the geometry of 
a charged black hole [5, 9] 



ds' 



with 



gMNdx M dx N 



f = l 



„2d-2 



M 



A, 



R 2 dr 2 
^) + 72-— ( 2 ) 



„rf— 2 



r d-2 



(3) 



ro is the horizon radius determined by the largest positive 
root of the redshift factor 



/M = o, 



M = rS 



„d-2 



(4) 



The caveat is there may exist certain operator dimensions or 
charges which do not arise in a consistent gravity theory with 
UV completion. 

The behavior at exactly zero frequency Gr(u) = 0, k) is con- 
trolled by UV physics, not by the IR CFT. 



7 It is defined so that for a typical supergravity Lagrangian it is a 
constant of order O(l). gF is related in the boundary theory to 
the normalization of the two point function of J M . 
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and 



B. AdS2 and scaling limits 



_ 9fQ 2{d-2) 



The geometry (2)-(3) describes the boundary theory 
at a finite density with the charge, energy and entropy 
densities respectively given by 



P = 



2(d-2) Q 



c d K 2 R d ^g F 7 



d-l M 
~2^~B J + T ' 



2tt (r \ d - 1 



n 2 \R 



(6) 
(7) 



The temperature of system can be identified with the 
Hawking temperature of the black hole, which is 



T = 



dr 
4irR 2 



1 - 



(d-2)Q< 
dr 2d ~ 2 



(8) 



In this paper we will be mostly interested in the be- 
havior of the system at zero temperature, in which limit 
the inequalities in (11) are saturated 



T = 



r = r, and M = ■ (16) 



Note that the horizon area remains nonzero at zero tem- 
perature and thus this finite charge density system has a 
nonzero "ground state" entropy density 8 , which can be 
expressed in terms of charge density as 



s = (27re d ) p 



(17) 



In the zero temperature limit (16) the redshift factor 
/ in (2) develops a double zero at the horizon 



/ = d(d _i) (!!_£=)! 



(18) 



and in (5) corresponds to the chemical potential. It can 
be readily checked from the above equations that the first 
law of thermodynamics is satisfied 



de = Tds + [idp 



(9) 



Note that Q has dimension of [L] d 1 and it is conve- 
nient to parameterize it as 



Q 



d-2 



(10) 



by introducing a length scale r*. In order for the met- 
ric (2) not to have a naked singularity one needs 



2{d - 1 rl 

~ d-2 * 



r > r* 



fin 



In terms of , the expressions for charge density p, chem- 
ical potential /x, and temperature can be simplified as 



P ~^\RJ e d ' 



/' 



d(d — 1) r* / r 
d-2 R 2 



ro 



d-2 



T=-^Ul 



AttR 2 

where we have introduced 

e d = 



2d-2 



2d-2 



9f 



v / 2d{d- 1) 



(12) 
(13) 
(14) 

(15) 



Note that r* can be considered as fixed by the charge 
density of the boundary theory. 



As a result, very close to the horizon the metric becomes 
AdS2 x R d_1 with the curvature radius of AdS2 given by 



R 2 



1 



-.R 



(19) 



y/d(d-l)' 

More explicitly, considering the scaling limit 
r-r* = A-^, t = \~ 1 T, A -> with C, r finite (20) 



we find that the metric (2) becomes a black hole in 

2 



AdSn x R d_1 



with 



(21) 



(22) 



The scaling limit (20) can also be generalized to finite 
temperature by writing in addition to (20) 



R 2 

ro — r* = A — - with (o finite 
Co 



(23) 



after which the metric becomes 



ds 2 = ^ 
Ub ^ 2 



1 - J) dr 2 



dC 



1 < 



f 2 dx 2 (24) 



Since the semiclassical gravity expression for the entropy is valid 
in the large N limit, one only needs "ground state degeneracy" 
in the N — > oo limit. 
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with 



and a temperature (with respect to r) 

1 



T 



2^Co 



(25) 



(26) 



Note that in the scaling limit (20), finite r corre- 
sponds to the long time limit of the original time coordi- 
nate. Thus in the language of the boundary theory (21) 
and (24) should apply to the low frequency limit 



1 £-0 



(27) 



where u> is the frequency conjugate to t. 



C. Emergent IR CFT 

One expects that gravity in the near-horizon AdS2 re- 
gion (21) of an extremal charged AdS black hole should 
be described by a CFTi dual. Little is known about 
this AdS2/CFTi duality 9 . For example, it is not clear 
whether the dual theory is a conformal quantum me- 
chanics or a chiral sector of a (1 + l)-dimensional CFT. 
It has been argued in [6] that the asymptotic symmetry 
group of the near horizon AdS2 region is generated by 
a single copy of Virasoro algebra with a nontrivial cen- 
tral charge, suggesting a possible description in terms of 
some chiral 2d CFT 10 . Some of the problems associated 
with AdS2, such as the fragmentation instability and the 
impossibility of adding finite-energy excitations [12] are 
ameliorated by the infinite volume of the factor in 

the geometry (21). 

The scaling picture of the last subsection suggests that 
in the low frequency limit, the d-dimensional boundary 
theory at finite charge density is described by this CFTi, 
to which we will refer below as the IR CFT of the bound- 
ary theory. It is important to emphasize that the confor- 
mal symmetry of this IR CFT is not related to the mi- 
croscopic conformal invariance of the higher dimensional 
theory (the UV theory) which is broken by finite charge 
density. It apparently emerges as a consequence of col- 
lective behavior of a large number of degrees of freedom. 

In section III we will elucidate the role of this IR CFT 
by examining the low frequency limit of two-point func- 
tions of the full theory. Our discussion will not depend 



9 For previous work on AdS2/CFTi correspondence from other 

decoupling limits, see e.g. [11] and its citations. 
10 The central charge is proportional to the volume of the d — 1- 
dimcnsional transverse space and is thus infinite for (21). To 
have a finite central charge one could replace R d_1 in (21) by a 
large torus. 



on the specific nature of the IR CFT, but only on its 
existence. In Appendix D we give correlation functions 
for a charged scalar and spinor in the IR CFT as cal- 
culated from the standard AdS/CFT procedure in AdS2 
[19]. They will play an important role in our discussion 
of section III. 



III. LOW FREQUENCY LIMIT OF RETARDED 
FUNCTIONS 



In this section we elucidate the role of the IR CFT 
by examining the low frequency limit of correlation func- 
tions in the full theory. We will consider two-point re- 
tarded functions for simplicity leaving the generalization 
to multiple-point functions for future work. We will 
mostly focus on zero temperature. 

Our discussion below should apply to generic fields in 
AdS including scalars, spinors and tensors. We will use a 
charged scalar for illustration. The results for spinors will 
be mentioned at the end with calculation details given 
in Appendix A. Vector fields and stress tensor will be 
considered elsewhere. 

Consider a scalar field in AdS^+i of charge q and mass 
to, which is dual to an operator O in the boundary CFT^ 
of charge q and dimension 



A 



2 i? 2 



d 2 



(28) 



In the black hole geometry (2), the quadratic action for 
(j) can be written as 

S = - J d d+1 xy/^j [(D M (f>)*D M (f) + to 2 0*0] (29) 

with 

D M (f> = (d M - iqA M )<t> . (30) 

Note that the action (29) depends on q only through 

fi q = fj,q (31) 

which is the effective chemical potential for a field of 
charge q. Writing 11 



d d k 



(2tt 



(32) 

the equation of motion for 0(r, k^) is given by (below 
k 2 = \k\ 2 ) 



1 



--d r {^g rr dr4>) + (g"(k 2 - u 2 ) + m 2 ) = 



(33) 



For simplicity of notation, we will distinguish <fr(r, x 11 ) from its 
Fourier transform <f>(r, by its argument only. 
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where various metric components are given in (2) and i^t _> 0. Note that since the definition of £ involves u, 

the matching will reshuffle the perturbation series in two 
regions. 

-gtt \ ' \ r d ~ 2 J J While the above scaling limit is defined for small u> > 

0, all our later manipulations and final results can be 
analytically continued to generic complex values of |cj| <C 
1. 



u{r) 



W + pL q 1 



d-2 



(34) 



In (3) we have chosen the gauge so that the scalar 
potential is zero at the horizon. As a result At — * fJ> for 
v — ► oo and u(r — » oo) — > + This implies that 
to should correspond to the difference of the boundary 
theory frequency from fi q . Thus the low frequency limit 
really means very close to the effective chemical potential 
/V 

The retarded Green function for O in the boundary 
theory can be obtained by finding a solution <f> which 
satisfies the in-falling boundary condition at the horizon, 
expanding it near the boundary as 



0, 



„A-d 



B(k M ): 



-A 



and then [10] 



G R (k^) = K 



B{K) 



(35) 



(36) 



where K is a positive constant which depends on the 
overall normalization of the action, and is independent 
of k u . 



1. Inner region: scalar fields in AdS2 

The scaling limit (37), (39) is in fact identical to that 
introduced in (20) (with u) replacing A) in which the met- 
ric reduces to that of AdS2 xR d_1 with a constant electric 
field. It can then be readily checked that in the inner re- 
gion at leading order, equation (33) (i.e. the equation 

for (jy'p) reduces to equation (D3) in Appendix D for a 
charged scalar field in AdS2 with an effective AdS2 mass 



(42) 



A single scalar field <j> in AdS^+i gives rise a tower of fields 
in AdS2 labeled by the spatial momentum k. From 
the discussion of Appendix D, the conformal dimension 
for the operator in the IR CFT dual to 4>j: is given by 



A. Low frequency limit 



At T = 0, expanding (36) in small to is not straight- 
forward, as the lo — > limit of equation (33) is singu- 
lar. This is because g tf has a double pole at the horizon. 
As a result, the w-dependent terms in equation (33) al- 
ways dominates sufficiently close to the horizon and thus 
cannot be treated as small perturbations no matter how 
small cu is. To deal with this we divide the r axis into 
two regions 



Inner: 
Outer: 



< r — r.. 



for e < C < oo (37) 
(38) 



and consider the limit 

lo — ► 0, ( = finite, e 



0. 







(39) 



Using £ as the variable for the inner region and r as that 
for the outer region, small w perturbations in each region 
can now be treated straightforwardly, with 



inner 



outer 



<MC) = ^ 0) (0+^ 1) (0 + --- (40) 



(!)/ 



.(41) 



We obtain the full solution by matching <fij and <fro in 
the overlapping region, which is — > with r — r* = 



Note that momentum conservation in M. d 1 implies that 
operators corresponding to different momenta do not 



0t(i)0£,(O)\ x5(k-k')t 



-2S k 



(44) 



To compute the retarded function (36) for the full the- 
ory, we impose the boundary condition that q^-p should 
be in-falling at the horizon. Near the boundary of the 



inner region (AdS2 region), i.e. Q 
then be expanded as (see (D5)) 12 



UJR7, 



0, 



r — r* 



(i + o(0) 

(1 + O(0) (45) 



R 2 X 2+^ 



The coefficient Gk{^) of the second term in (A32) is pre- 
cisely the retarded Green function for operator in the 



For convenience for matching to outer region we have taken a 



,(()) 



below. The calculation 



specific choice of normalization for Vj 
of retarded function (36) does not depend on the choice of nor- 
malization. 
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IR CFT. From (D10) it can be written as 
Gk (w 



T(2v k )T (i - v k - iqe d ) 

(46) 

with v k given by (43). Equation (45) will be matched to 
the outer solution next. 



amounts to choosing a specific special solution of the 
homogeneous equation. Note that it is important that 
the equations are linear. Given that higher order terms 
in (49) are uniquely determined by rj±^ , to match the full 
solution (f>o to the inner region it is enough to match the 
leading order term which we have already done. We thus 
conclude that perturbatively 



2. Outer region and matching 



4>0 = V+ + Gk(u)r). 



(50) 



The leading order equation in the outer region is ob- 
tained by setting lu = in (33). Examining the resulting 
equation near r — ► r* , one finds that it is identical to the 
inner region equation for <pf in the limit £ — > 0. It is 
thus convenient to choose the two linearly-independent 
solutions rj±^ in the outer region using the two linearly 

independent terms in (45), i.e. rj^ are specified by the 
boundary condition 



i?2 



. (47) 



The matching to the inner region solution (45) then be- 
comes trivial and the leading outer region solution epQ 
can be written as 



(48) 



with Q k [uj) given by (46). 

One can easily generalize (48) to higher orders in lo. 
The two linearly independent solutions to the full outer 
region equation can be expanded as 



(o) 



(i) 



2 (2) 



(49) 



where higher order terms r)± , n > 1 can be obtained 
using the standard perturbation theory and are uniquely 
specified by requiring that when expanded near r = r*, 
they do not contain any terms proportional to the zeroth 
order solutions. Each of the higher order terms satis- 
fies an inhomogenous linear equation. The requirement 



3. Small lu expansion of Gr 

We first look at the retarded function at lu = 0. At 
lu = the inner region does not exist and the outer re- 
gion equation reduces to that satisfied by (j^Q ■ In (48) 
we have chosen the normalization so that at to 

„(o) 



0, 

?7 < ^ , ' ) . For real u k , this follows from the fact that 



r) + gives the regular solution at 



When Vk is pure 



imaginary (i.e. when q is sufficiently large) we will define 
the branch of the square root by taking m 2 — > m 2 — ie 
so that Vk = —i\k with Afe positive. Then rj^ is the 
in-falling solution at the horizon as is required by the 
prescription for calculating retarded functions. Expand- 
(0)/ \ 

mg rj\ (r) near r — > oo as 



v f>(r, k) = af\k)r A - d (1 + • • •) + b^ik)^ (1 + • • •) 

(51) 

then from (36) we find that 



G r {lo - 0, k) = K - 



h (0) 

w 



(52) 



Now consider a small nonzero to. Expanding various 
functions in (49) (n > 1) near r — > oo as 

^(r, k) = a£\k)r*-* (1 + • • •) + b { ]!\k)r- A (1 + • • •) , 

(53) 

from (50) and (36) we find that for small co 



Gr(lo, k) = K 



bf + Ujb^ + 0(LU 2 ) + g k {L0) (b^ + Lub W + 0(U 2 )) 

+ ua ( l ] + 0{uj 2 ) + g k (uj) (a^ + + 0(w 2 )) 



(54) 



Equation (54) is our central technical result. In next few 
sections we explore its implications for the low energy 
behavior of the finite density boundary system. While 
its expression is somewhat formal, depending on various 
unknown functions a±\k), b±\k) which can only be ob- 



tained by solving the full outer region equations order by 
order (numerically), we will see that a great deal about 
the low energy behavior of the system can be extracted 
from it without knowing those functions explicitly. 
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B. Generalization to fermions 

Our discussion above only hinges on the fact that in 
the low frequency limit the inner region wave equation 
becomes that in AdS2- It applies also to spinors and 
other tensor fields even though the equations involved are 
more complicated. In Appendix A we discuss equations 
and matching for a spinor in detail. After diagonalizing 
the spinor equations one finds that eigenvalues of the re- 
tarded spinor Green function (which is now a matrix) are 
also given exactly by equation (54) with now Q k (uj) given 
by equation (A34), which we copy here for convenience 



T(-2v k )T(± + v k -iqe d ) 
Y{2v k ) r (1 

ikR \ R 



iqe d - Vk 



ikR 



with 



R 2 - iqe d + v k 



k 2 R 2 



"k - —pr 



{2uY Uk (55) 



(56) 



Note the above scaling exponent can also be expressed 
as (using (13)) 



Vk 



9FQ 



2m 2 R 2 d(d-l)k 2 



V2d(d-1) V 9W (rf-2) 2 A 



(57) 



The conformal dimension of the operator in the IR 
CFT is again given by 



• For both scalars and spinors, e' 7fc (and thus c(k)) 
always lies in the upper-half complex plane. 

• For scalars e ZJk+2m ' /k always lies in the lower-half 
complex plane, while for spinors e %lk+2lx%Vk always 
lies in the upper-half complex plane. 



For v k € (0, i), 



spinor 
scalar : 



7T - 7^ > 2-KV k 
7T - 7fc < 2-KV k 



(61) 



For pure imaginary v k = —i\ k (X k > 0), the ratios 
in (D25) and (D26) of Appendix D become real and give 



-2lT\ k 



,-2ivqe d 



Cfe 



and 



e 27rAfe _ e -2irqe d 



e -2nX k 1 e -2Trqe d 



< e 



-47rA fc 



spmor 



(62) 



Cfe 



o 2 7T A f, 



-2irqe d 



> e" 47rAfc scalar . (63) 



It is also manifest from the above expressions that 
|c(fc)| 2 < 1 for both scalars and spinors. 

Also note that for generic v k , Q k {uj) and accordingly 
Gr(uj, k) in (54) have a logarithmic branch point at u> = 

0. We will define the physical sheet to be 6 £ (— f , 

1. e. we place the branch cut along the negative imaginary 
axis. This choice is not arbitrary. As discussed in Ap- 
pendix D, when going to finite temperature, the branch 
cut resolves into a line of poles along the negative imag- 
inary axis. 



1 

Ok = 2 + ■ 



(58) D. Renormalization group interpretation of the 

matching 



C. Analytic properties of Q k 

The analytic properties of Gk(&) wm pl&y au important 
role in our discussion of the next few sections. We collect 
some of them here for future reference. Readers should 
feel free to skip this subsection for now and refer back to 
it later. 

We first introduce some notations, writing 

g k (u) = c(k)u? Vk , c(k) ee |c(fc)|e^ fc (59) 

where c(k) denotes the prefactor in (55) for spinor and 
that in (46) for scalars. 

For real v k , the ratios in (D25) and (D26) of Ap- 
pendix D become a pure phase and we find that 

!arg (T(-2v k ) (e" 27rWfc - e" 27 ^)) spinor 
arg (r(-2i/ fc ) ( e ~ 2m » k + e - 2 ^)) scalar 

(60) 

It can be readily checked by drawing e 27rw *= and e 27TC i ed 
on the complex plane that the following are true: 



The matching procedure described above has a natural 
interpretation in terms of the renormalization group flow 
of the boundary theory. The outer region can be inter- 
preted as corresponding to UV physics while the inner 
AdS2 region describes the IR fixed point. The matching 
between in the inner and outer regions can be interpreted 
as matching of the IR and UV physics at an intermediate 
scale. More explicitly, coefficients , from solving 
the equations in the outer region thus encode the UV 
physics, but Gk{u) is controlled by the IR CFT. 

In this context u can be considered as a control pa- 
rameter away from the IR fixed point. Equation (54) 
then shows a competition between analytic power cor- 
rections (in to) away from the fixed point and contribu- 
tion from operator 0g. In particular when v k > \ {i.e. 
5 k > 1), 0£ becomes irrelevant in the IR CFT and its 
contribution becomes sublcading compared to analytic 
corrections. Nevertheless the leading non- analytic contri- 
bution is still given by Gk{w) and as we will see below in 
various circumstances Q k does control the leading behav- 
ior of the spectral function and other important physical 
quantities like the width of a quasi-particle. 
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It is interesting to note the similarity of our match- 
ing discussion to those used in various black hole/brane 
emission and absorption calculations (sec e.g. [13]) 
which were important precursors to the discovery of 
AdS/CFT. The important difference here is that in this 
asymptotically- AdS case we can interpret the whole pro- 
cess (including the outer region) in terms of the dual field 
theory 



V. EMERGENT QUANTUM CRITICAL 
BEHAVIOR II: LOG-PERIODICITY 

In this section we examine the implication of (54) when 
the fk becomes pure imaginary. We recover the log- 
oscillatory behavior for spinors first found numerically 
in [3]. 



IV. EMERGENT QUANTUM CRITICAL 
BEHAVIOR I: SCALING OF SPECTRAL 
FUNCTIONS 

In this and the next two sections we explore the im- 
plications of equation (54) for the low energy behavior 
of the finite density boundary system. In this section we 
look at the behavior of (54) at a generic momentum for 

which i>k is real and a+ ' (fc) is nonzero. Imaginary 
will be discussed in section V and what happens when 
or? (k) = will be discussed in section VI. 

When Vk is real, the boundary condition (47) is real. 



(0) 



Since the differential equation satisfied by <^jq 

real, one concludes that both b ® and a? are re; 
implies that 

Im Gr(uj = 0, k) = 0, for real v k . 



is also 



(64) 



Similarly we can conclude that all coefficients in (53) are 
also real. Thus the only complex quantity in (54) is the 
Green function of the IR CFT, g k (cu). When a^ 0) (fc) is 
nonzero we can expand the denominator of (54) and the 
spectral function for O can be written at small u> as 

lmG R (uj,k) = G R {k,u) = 0)d o Im£ fc (w) + • • • oc u 2Vk . 

(65) 

with 



d 



bf 



(66) 



We thus see that the spectral function of the full theory 
has a nontrivial scaling behavior at low frequency with 
the scaling exponent given by the conformal dimension of 
operator 0£ in the IR CFT. Note that the fc-dependent 

prefactor in (65) depends on a±\b±^ and thus the met- 
ric of the outer region. This is consistent with the RG 
picture we described at the end of last section; the scal- 
ing exponent of the spectral function is universal, while 
the amplitude does depend on UV physics and is non- 
universal. By "universal" here, we mean the following. 
We can imaging modifying the metric in the outer region 
without affecting the near horizon AdS2 region. Then 



t , u ± will change, but the exponent v k will remain 



the same. The real part of Gr is dominated by a term 
linear in ui when v k > \ and is non-universal, but the 
leading nonanalytic term is again controlled by Qk(w). 



A. Log-periodic behavior: complex conformal 
dimensions 

When the charge q of the field is sufficiently large (or 



m too small) 



f 



scalar : m 2 R\ + - < q 2 e\ 
spinor : 



(67) 



m R 2 <q e d 
there exists a range of momenta 



k 1 < ki = 



for which v k is pure imaginary 



rl 




w 












W 





scalar 



spinor 



(68) 



Ve3 



m 2 R 2 



m 2 R 2 



scalar 
spinor 



(69) 

We have chosen the branch of the square root of v k by 



taking m 



it. The effective dimension of the op- 



erator 0£ in the IR CFT is thus complex. Following [3] 
we will call this region of momentum space (68) the os- 
cillatory region. 13 For spinors we always have m 2 > 
and the existence of the oscillatory region requires q ^ 0. 
For scalars, equation (67) can be satisfied for q = for 
m 2 in the range 



d 2 n2 d(d - 
4 4 



(70) 



where the lower limit comes from the BF bound in 
AdSd+i and the upper limit is the BF bound for the 
near horizon AdS2 region. 

For a charged field an imaginary v k reflects the fact 
that in the constant electric field (22) of the AdS2 region, 
particles with sufficiently large charge can be pair pro- 
duced. It can be checked that equations (67) indeed coin- 
cide with the threshold for pair production in AdS2 [14]. 

With an imaginary the boundary condition (47) for 

?7± is now complex. As a result 6+ ' / a+ is complex and 



ImG R (u = 0,fc) ^ 



(71) 



13 Note that the oscillatory region appears to be different from the 
Fermi ball discussed in [2]. 
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Thus there are gapless excitations (since ui = 0) for a 
range of momenta k < k . This should be contrasted 
with discussion around (64). 

The leading small u> behavior (54) is now given by 



bf + b { ° ] ] c(k)m- 2tX « 
<4 0) +a ( ° ) c(k)u- 2iX " 



+ 0{lo) (72) 



where c{k) was introduced in (59). Note that here 



l(0) 



(b 



(oh 



(73) 



since r/^ = (ri { " > )* at the horizon and that the differ- 
ential equation the r\ satisfy is real. Equation (72) is 
periodic in logw with a period given by 



(0) 



Tk 



7T 



(74) 



In other words (72) is invariant under a discrete scale 
transformation 



'CO, 







(75) 



We again stress that while the retarded function (and 
the spectral function) depends on UV physics (i.e. solu- 
tions of the outer region) , the leading nonanalytic behav- 
ior in ui and in particular the period (74) only depends on 
the (complex) dimension of the operator in the IR CFT. 

Here we find that the existence of log-periodic behavior 
at small frequency is strongly correlated with (71), i.e. 
existence of gapless excitations. It would be desirable 
to have a better understanding of this phenomenon from 
the boundary theory side. 



B. (In)stabilities and statistics 

It is natural to wonder whether the complex expo- 
nent (69) implies some instability. We will show now 
that it does for scalars but not for spinors. The scalar 
instability arises because the scalar becomes tachyonic 
in the AdS2 region due to the electric field or reduced 
curvature radius. At zero momentum this is precisely 
the superconducting instability discussed before in [15- 
18]. 14 That the log-oscillatory behavior does not imply 
an instability for spinors was observed before in [3] by 
numerically showing there are no singularities in the up- 
per half w-plane. Below we will give a unified treatment 
of both scalars and spinors, showing that the difference 



14 It was noted before in [17, 18] that there could be an instability 
even for neutral scalar in the mass range (70). The condensation 
of such a neutral scalar field does not break the U(l) symmetry 
and thus is distinct from the superconducting instability of a 
charged field. It would be very interesting to understand better 
the boundary theory interpretation of such an instability. 



between them can be solely attributed to statistics even 
though we have been studying classical equations. 

The spectral function following from (72) can be writ- 
ten as 



lmG R (w,k) w > l-|c(fc)| : 



ImG fl (w = 0,fc) 



|c(fc)|e 



•X\2 



co 1 - \c(k)\ 2 e 4 ^ 



(76) 



|c(fc)|e'" 



v2n\ k p iX 12 



where we have introduced 



X = lk -2a-2\ k logM 



,(0) 



i«?V 



(77) 



and 7^ was defined in (59). In the boundary theory re- 
tarded Green functions for bosons are defined by com- 
mutators, while those for fermions by anti-commutators, 
which implies that for u> > 

scalars: Im G r (-lj, k) < 0, Im G r (uj, k) > (78) 
spinors: Im G r (~lu, k) > 0, Im G r (lu, k) > (79) 

Applying (78) and (79) to (76) requires that 

scalars : |c(fc)| 2 < 1, |c(fc)|e 27rAfc > 1 (80) 
spinors : |c(fc)| 2 < 1, |c(fc)|e 27rAfc < 1 (81) 

which arc indeed satisfied by (62)-(63). It is important to 
stress that in the bulk we are dealing with classical equa- 
tions of scalars and spinors and have not imposed any 
statistics. However, the self-consistency of AdS/CFT im- 
plies that classical equations for bulk scalars and spinors 
should encode statistics of the boundary theory. 

We now examine the poles of (72) in the complex lu- 
plane, which is given by 



1 



|c(fe)|e 



2\ k iX 



0, with lu — 



JO 



(82) 



(82) implies a series of poles located along a straight 
line with angle 8 C (in the expression below the integer 
n should be large enough for our small u approximation 
to be valid) 



7t -2c-(2n+l), 



-,i0 r 



2X k 



Equations (80) and (81) then imply that 

scalars : 6 C € (0, tt) 
spinors : 9 C > tt . 



-J-log|c(A;)| ■ 
(83) 

(84) 
(85) 



Thus the poles for scalars lie in the upper half w-planc 
[18] while those for spinors are in the lower half w-planc. 
Poles on the upper half w-plane of a retarded Green func- 
tion on the one hand implies causality violation. On the 
other hand from equation (36) it implies that there exist 
normalizable modes which have a frequency with posi- 
tive imaginary part. This leads to a mode exponentially 
growing with time and thus an instability in the charged 
black hole geometry. See figure 1 for illustration of the 
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locations of poles for scalars and spinors and their move- 
ment as k is varied in the range (68). 

The instability for a scalar can also be understood in 
terms of classical superradiance. To see this let us go 
back to equation (45) which for i>\. = —i\k the first term 
can be interpreted as an incident wave into the AdS2 
region with the second term the reflected wave. Thus 
the reflection probability is given by 



|c(fc)| 5 
\c(k)\ 



LO > 

lu < 



(86) 



Equation (80) then implies |<?fc(u>)| > 1 for co < 0. Re- 
call that our lu is defined to be the deviation from the 
effective chemical fi q = qfi. Thus lu < agrees with the 
standard frequency region for superradiance. As men- 



tioned earlier, an imaginary v corresponds to the parame- 
ter regime where charged particles can be pair-produced. 
While both scalars and spinors are pair-produced, the 
superradiance of scalars can enhance the pair production 
into a classical instability. The produced scalar particles 
are trapped by the AdS gravitational potential well (of 
the full geometry), and return to the black hole to induce 
further particle production. In contrast, since the reflec- 
tion probability for a spinor falling into the black hole is 
smaller than 1, after a few bounces back and forth the 
paired produced spinor particles should fall back into the 
black hole. 

In our context, the fact that a scalar superradiates 
while a spinor does not (i.e. equations (80)-(81)) can 
be seen as a consequence of statistics of the operator in 
the boundary theory. 




FIG. 1: The motion of poles of the Green functions (72) of spinors (left) and scalars (right) in the complex frequency plane. 
For illustration purposes we have chosen parameters and rescaled |o;|( — > |u;| # with small #) to give a better global picture. 
The poles are exponentially spaced along a straight line (dotted line) with angle 8 C given by (83). There are infinitely many 
poles, only a few of which are shown. Left plot: The black dashed line crossing the origin corresponds to the value of 8 C at 
k = k (see (68)): the boundary of the oscillatory region. As k — > k , most of the poles approach the branch point lj = except 
for a finite number of them which become quasi-particle poles for the Fermi surfaces at larger values of k. The color dashed 
lines in the right half indicate the motion of poles on another sheet of the complex frequency plane at smaller values of k (see 
the end of sec. IIIC for the choice of branch cut in the w-plane). Right plot: the two dashed lines correspond to k — (upper 
one) and k = k (lower one). Again most of the poles approach the branch point u = as k — > k . These plots are only to be 
trusted near lu = 0. 



VI. EMERGENT QUANTUM CRITICAL 
BEHAVIOR III: FERMI SURFACES 



A. Quasi-particle-like poles 



We now consider what happens to (54) when a+\k) 
in the expansion of (51) is zero. This can only occur for 
real at discrete values of k, at which values the wave 
function rj^ becomes normalizable. The possible exis- 
tence of such fc's can be visualized at a heuristic level by 
rewriting (33) with lu = in the form of a Schrodingcr 
equation, and noticing that for a certain range of momen- 
tum the Schrodingcr potential develops a well which may 
allow normalizable "bound states" . See Appendix B for 



details and also the similar story for spinors. For which 
values of k such bound states indeed occur can then be 
determined by solving (33) (with lu = 0) numerically. 



Suppose that 
from (52), for k > 



G R (k,u = 0) 



(0) 



has a zero at k = ki 



d k af{k F ) 



k i = k — kf 



Then 



(87) 



Since v k is real, a + , b , are all real (see e.g. discussion 
around (64)). Thus lmGn(uj = 0, k) is identically zero 
around k F , but the real part ReGfl(w = 0, k) develops a 
pole at k — k F . Now turning on a small lo near k F , we 
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then have to leading order 



d k a^ (k F )kj_ + ujo}+> (k F ) + a K "> {k F )Q kF (, 



(0)/ 



hi 



(88) 



v F 



aL 0) (^) 



|c(MI 



^(M 



(89) 



The term linear in u> in the downstairs of (88) can be 
omitted if v kF < i . For v kp > | we should still keep the 
term proportional to u 2l,k F , since it makes the leading 
contribution to the imaginary part. Note that the quan- 
tities in (89) are all real and as we will discuss below they 
are all positive. 

The coefficients v F , hx,2 may be expressed in terms of 

integrals of bound state wave function 77+ at k = k F by 
perturbing the Schrodinger problem at to = 0, k = k F 
in cj, k, similar to the demonstration of the Feynman- 
Hellmann theorem. We present the details of this analy- 
sis in Appendix C. In particular, from a combination of 
analytic and numerical analysis, we show that (for q > 0): 

1. For v kp > i, v F > for both scalars and spinors 
and in particular for a spinor v F < 1. 

2. For all v kp and for both scalars and spinors 

h 1 ,h 2 >0. (90) 

For scalars the above equation is established an- 
alytically in Appendix C2. For spinors it can be 
checked numerically, see fig. 8. 

Equation (88) leads to a pole in the complex-w plane 15 
located at 



uj c (k) = oj*{k) — iT(k) 



ik. 



•'I'F 



v F k ± - foe*"* {v F k±) 2vk *- v kF > \ 



with residue at the pole given by 



Notice that both the real and imaginary part of the pole 
go to zero as k± — > 0. Thus (91) leads to a sharp quasi- 
particlc peak in the spectral function IuiGr(uj, k) in the 
limit fcx — > with a dispersion relation 



where in the second line we have used (59) and introduced an( j 



u*(k) oc k\ with z 



r(fc) cx k{ with 5 




v kp < 
v kF > 




v kp < 
v kp > 



(93) 



(94) 



Note that when v kp < i, the pole follows a straight 
line as k± is varied. More explicitly, 



arg(w c ) 



2^ 

1k f 



k ± >0 
k± < 



(95) 



and the width T is always comparable to the frequency 

, , 16 



U!*(k) 



tan r 



r(fe) 



In contrast, for v kp > i, - ^ 



const 



as k\ 



(96) 



0. See 



fig. 3 for examples of the motion of a spinor pole in the 
complex w-plane. 

At v kp = i, a^ 1 ' and G k {uj) in equation (88) are diver- 
gent. Both of them have a simple pole at v kp = |. The 
pole in Q k can be seen explicitly from the factor r(— 2v k ) 

in (46) and (55). The pole in can be seen from the 
discussion (C8)-(C9) in Appendix C. The two poles can- 
cel each other and leave behind a finite to log u> term with 
a real coefficient. More explicitly, we have 



Gt 



hi 



k± + c\oj log u) + c\lo 



(97) 



where c\ is real and ci is complex. 

Similar logarithmic terms appear for any v kp 
t|, n G Z_|_. For example, at v kp = 1, one finds that 



fc_L - -^UJ + c 2 ijJ 2 logcj + C 2 UJ 2 
with £2 real and c 2 complex. 

B. A new instability for bosons 



(98) 



Z = 



u c hi 
2u kp k± 

-hiv F 



oc k 



v kF < 
v kp > 



The above discussion applies identically to both scalars 

(92) 

and spinors with their respective parameters. We now 



Recall the discussion at the end of sec. Ill C for the choice of 
branch cut in the w-planc. 



Note the concept of a width is only operationally meaningful 
when the pole lies in the physical sheet. See the end of sec. Ill C 
for our choice of the physical sheet. 
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show that for spinors the pole (91) never appears in the 
upper half plane of the physical sheet, while for scalars 
it always lies in the upper half plane for k± < 0. The 
difference can again be attributed to the statistics of the 
corresponding boundary operators. 

First note that Im Gr obtained from (88) should again 
satisfy (78) and (79) (which follow from the statistics of 
the full boundary theory), leading to (using also (90)) 

scalars : sin7fc F > 0, sin (7^^ + 27r^ F ) < (99) 
spinors : s\n^k F > 0, sin (7^ + 2-Ku kp ) > (100) 

which indeed follow from discussion below (60). Equa- 
tions (99) and (100) are also consequences of (for u > 0) 

scalars: Im<? fc (-o;) < 0, hnQ k {iJ) > (101) 
spinors: JmQ k (-u) > 0, Imft(w) > (102) 

which follow from the Bose and Fermi statistics of the IR 
CFT. This gives a self-consistency check of the statistics 
of the full theory and its IR CFT. 

For z/fc F > i, applying equations (99)-(100) to the last 
line of (91) we find that for spinors the pole always lies in 
the lower half plane while for scalars the pole is always 
on upper half complex w-planc for fcj_ < 0. By using 
in addition (61) we again reach the conclusion that for 



spinors the pole never appears in the upper half plane 
of the physical sheet, while for scalar it always does for 
k± < 0. In fig. 2 we give a geometric picture to illustrate 
this. In fig. 3 we illustrate the motion of poles as k±_ is 
varied for spinors. 

The pole in the upper half plane for scalars when 
k± < is intriguing. As mentioned earlier poles in the 
upper-half cj-planc imply the existence of exponentially 
growing (in time) normalizable mode in AdS and lead to 
instability. Again as in the case of the oscillatory region 
the instability for scalars can be attributed to the Bose 
statistics of the boundary theory. This instability is cu- 
rious as it occurs for real v k and is thus distinct from 
the instability discussed in last section, which is associ- 
ated with an imaginary v k . In particular, it appears that 
this instability can exist in a parameter range where pre- 
viously considered superconducting instability does not 
occur. To understand the physical interpretation of such 
an instability we should examine the motion of the pole 
as k is decreased, in particular, whether it persists to 
k = 0. If it does, it seems likely the pole will have finite 
real part at k = 0, i.e. the growing mode also oscillates in 
time. Such an instability appears to be novel and we will 
leave its interpretation and a detailed study for future 
work. 
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FIG. 2: A geometric illustration that poles of the spinor Green function never appear in the upper-half oj-plane of the physical 
sheet, for two choices of Uk F < \. Depicted here is the io 2uk F covering space on which the Green function (88), with the uj/vf 
term neglected, is single- valued. The shaded region is the image of the upper-half o;-plane of the physical sheet. The pole lies 
on the line 1Vk F Q c = — 7fe F for k± > and on 2v^ F 9 c — n — jkp for k± < 0, which are indicated by the purple solid line in 
the figure. The triangle formed by dashed arrows and solid lines in the upper left quadrant gives the geometric illustration 
for the equation 7T — 7^ = axg(e 2n ™ k — e~ 27Tqed ) (following from the first equation of (60)), which makes it manifest that for 
k± < the pole lies outside the shaded region. In contrast for a scalar one needs to reverse the direction of the horizontal 
dashed line and the pole lies inside the shaded region. Similarly, the triangle in the lower right quadrant gives the illustration 
for — 7fc = arg(— e 2nl,Uk + e~ 27Tqed ) which is relevant for k± > 0. We also indicated the angles ot± which will be introduced and 
discussed in detail around (104). 



C. Fermi surfaces 

We now focus on spinors, for which equations (88) 
and (91) give analytic expressions for, and generalize to 



any mass m and charge q, numerical results of [3]. Ref- 
erence [3] focused on m = and a few values of charge 
q. There kp was found by studying the scaling behav- 
ior of quasi-particle peaks as they become sharper and 
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FIG. 3: Examples of the motion of the pole for a spinor as k± is varied (arrows indicating the directions of increasing k±). Left 
plot: Vk F < |, for which the pole moves in a straight line. The plot shows an example where the pole moves to another sheet 
of the Riemann plane for k± > (i.e. a + > -|); a± indicated there are introduced in (104). Right plot: Vk F > | for which the 
dispersion (real part of the pole) is linear. 



sharper as kp is approached. In our current discussion 
&f are found from "bound states" of the Dirac equation 
at co = (which needs to be found numerically). We 
found perfect agreement between two approaches. Plug- 
ging explicit values of kp into (88) and (91) also leads to 
almost perfect agreement with numerical plots of the re- 
tarded function in [3] including the scaling exponents (see 
ng. 4). 
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FIG. 4: Spinor Green's function G2 (as defined in Ap- 
pendix A) at k = 0.918 as a function of uj, computed nu- 
merically. We have chosen parameters r* = R = gp = q = 1, 
m = and d = 3 for which the Fermi momentum is kF — 
0.91853. The real and imaginary parts are shown in blue and 
orange dotted curves, respectively. Also shown is (88) (solid 
lines) with hi , hi computed numerically using the method of 
Appendix C and Vk given by (56). 

The sharp quasi-particlc peaks in the spectral function 
for spinors were interpreted in [3] as strong indications 
of underlying Fermi surfaces with Fermi momentum kp. 
The scaling behavior near a Fermi surface for the param- 
eters range considered there, which all have v kp < |, is 
different from that of the Landau Fermi liquid, suggesting 
an underlying non-Fermi liquid. Our current discussion 
allows us to obtain a "landscape" of non-Fermi liquids 
by scanning all possible values of m and q. Before do- 
ing that, let us first make some general comments on the 



scaling behavior (91)-(96) and the possible underlying 
non- Fermi liquids: 

1. The form of retarded Green function (88) and scal- 
ing behavior (91) again have a nice interpretation 
in terms of the RG picture described earlier; while 
the location of the Fermi momentum kp is gov- 
erned by the UV physics (just like in real solids), 
the scaling exponents (91) near a Fermi surface are 
controlled by the dimension of the corresponding 
operators £>£ (with |fc| = kp) in the IR CFT. Note 
that while the specific value of Vk F depends on kp 
as an input parameter, its functional dependence 
on kp is fixed by the IR CFT. 

2. For Vk F < \, the corresponding operator in the IR 
CFT is relevant 17 . From (96), the imaginary part 
of the pole is always comparable to the real part 
and thus the quasi-particle is never stable. Also 
note the ratio (96) depends only on the IR data. 
Another important feature of the pole is that its 
residue (92) goes to zero as the Fermi surface is ap- 
proached. In particular, the smaller v kF , the faster 
the residue approaches zero. 

3. When v kp > i, the corresponding operator in the 
IR CFT is irrelevant. The real part of the dis- 
persion relation (91) is now controlled by the an- 
alytic UV contribution and becomes linear with a 
"Fermi velocity" given by vp. Vf is controlled by 
UV physics and will have to be found numerically 
by solving the outer region equations (see fig. 7 and 
discussion of the next subsection). The imaginary 
part of the pole is still controlled by the dimension 
of the operator in the IR CFT. In this case the 
width becomes negligible compared with the real 
part as the Fermi surface is approached k± — > 0, 
so the corresponding quasi-particle becomes sta- 
ble. Furthermore the quasiparticle residue is now 



Recall that the dimension in the IR CFT is &k F = 2 F ^k F ■ 
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non- vanishing (92). Note the scaling exponent of 
the width T is generically different from the in- 
dependence of the Landau Fermi liquid. 

4. When = ^, the corresponding operator in the 
IR CFT is marginal. Now the retarded function is 
given by (97), in which case the imaginary part 
of the pole is still suppressed compared to the 
real part as the Fermi surface is approached, but 
the suppression is only logarithmic. The quasi- 
particlc residue now vanishes logarithmically as the 
Fermi surface is approached. Remarkably (97) 
is precisely of the form postulated in [8] for the 
"Marginal Fermi Liquid" to describe the optimally 
doped cuprates. In [8] the term "marginal" referred 
to the logarithmic vanishing of the quasi-particle 
weight (residue) approaching the Fermi momen- 
tum. Here we see that this term "marginal" is in- 
deed perfectly appropriate. 

5. At Vk F = 1, from (98) the retarded Green function 
resembles that of a Landau Fermi liquid with the 
real part of the quasi-particle pole linear in k± and 
the imaginary part quadratic in k±. Equation (98), 
however, has a logarithmic term aj 2 logw with a 
real coefficient which is not present for a Landau 
Fermi liquid 18 . In particular, this leads to a curious 
particle-hole asymmetry with a difference in width 
for a hole and a particle given by 19 



T(w* < 0) - r(w* > 0) = tt5 2 w 2 



(103) 



6. Notice from (93) that for all values of Vk F , z > 1. 
This is consistent with an inequality proposed by 
Scnthil in [7] for a critical Fermi surface. 20 In 
the notation of [7] our Green function (88) also 
has scaling exponent a = 1, with a defined as 
Gr{\ z u, AfcjJ = \- a G R {uj 1 k 1 _). Thus we also 
have z > a for all cases, consistent with the other 
inequality proposed in [7]. 

7. For Vk F < \, equation (88) (with the linear term in 
lu omitted) exhibits a particle-hole asymmetry. To 
characterize this it is convenient to define 



n + 



-0 c (fcj. > 0) = 

e c (k x < o) - 7T 



ik F 

1vk F 



(104) 



with 9 C introduced in (95). a+ gives the angular 
distance away from the positive real u-axis for the 
pole (i.e. particle- type excitations) at k± > 0, while 
a_ gives the angular distance away from the neg- 
ative real w-axis for the pole (i.e. hole- type excita- 
tions) at < 0. From fig. 2 and its caption, a± 
are always positive. Their values determine how 
close the corresponding pole stays to the real axis 
and hence indicate the sharpness of the resulting 
peak in the spectral function (i.e. imaginary part 
of Gr). In particular, when a± exceeds ^ the corre- 
sponding pole moves to the other Riemann sheet 21 . 
The particle-hole asymmetry is then reflected in the 
relative magnitudes of a±. From equation (60) one 



finds that 2vk F a.- = arg(l 



-2-xqe d — 2-Kiv k 



f ) which 



implies a_ quickly becomes small when q increases 
as can also be visualized from fig. 2 since the hori- 
zontal line of the upper triangle becomes short. In 
contrast, with some thought one can conclude from 
fig. 2 that a+ can only be small when Vk F is close to 
i, where a_ is also small. Away from Vk F = |, one 
will generically (except for small q) have a particle- 
hole symmetry with a sharp peak on the hole side 
(i.e. small a_), but a broad bump on the particle 
side (not so small a+). This was indeed what was 
observed in [3] numerically. The above qualitative 
features will be confirmed by explicit numerical cal- 
culations represented in fig. 6. 



8. Non-Fermi liquids have been described previously 
by coupling a Fermi surface to a propagating 
bosonic mode, such as a transverse magnetic exci- 
tation of a gauge field, e.g. [22-29]. The forms of 
the fermion Green's functions thus obtained all fit 
into the set of functions we have found in (88). An 
important difference, however, is that each of these 
analyses required a small parameter 22 to control 
the perturbation theory, and the range of frequen- 
cies (or temperatures) over which the non-Fermi 
liquid behavior is relevant is parametrically small 
in the control parameter. As a result the non-Fermi 
liquid behavior will only be visible at extremely 
low temperatures. In our non-perturbativc calcula- 
tion, this range is order unity, and so our non- Fermi 
liquids may be considered robust versions of these 
previously-identified phases. 



18 It is amusing to note that a Landau Fermi liquid in (2 + l)-d has 
a logarithmic correction of the form u) 2 log u>, but with a pure 



1!) 



imaginary coefficient (see e.g. [20]). 
The above asymmetry can also be expressed as 

p -2nqe d 



nw<o) 
r( W „>o) 



20 Senthil derived the inequality z > 1 by approaching a non-Fermi 
liquid at the critical point of a continuous metal-insulator transi- 
tion from the Landau Fermi liquid side and requiring the effective 
mass should not go to zero as the critical point is approached. 



9. The expression (88) for retarded function near a 
Fermi surface and the matching procedure from 
which (88) was derived suggest the following ef- 
fective action of the UV On and IR Oi part of an 



21 An example is given in the left plot of fig. 3. 

22 e.g. a ~ in [22], 1/N in [23], the parameter x in the gauge 
boson dispersion in [24]... 
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operator O (below \k F \ = kp) 

S = J dcjdk6 u Y,(uj,kx)Ou 
+ J dudk F D(uj,k F )6 u (uj,k F yO I (u,k F ) + h.c. (105) 

where S(w, fcj_) represents the kinetic term for Ojj 
and -D(o;, /cp) denotes coupling between Oi and 
Ojj. Both £ and D are controlled by UV physics. 
They are assumed to be real and depend analyti- 
cally on uj. We assume that the dynamics of the 
IR operator Oi(cu,k F ) is controlled by the IR CFT 
with two point function given by 

(106) 

The action then implies that after summing a geo- 
metric series, the full correlation function of Ojj is 
then given by 



G R (u,k A 



Z(uj 7 k ± )+D 2 {Lu,k F )g kF (uj) 



(107) 



which has the form of (88). Thus, action (105) gives 
a phenomenological model of the small-frequency 
matching procedure described in this paper as a 
coupling between UV and IR degrees of freedom. 
This argument is not dissimilar to that taken in [8] 
to obtain the "marginal Fermi liquid" . It might also 
be possible to reinterpret the discussion of [22-29] 
in this language. 



D. The zoo of non- Fermi liquids from gravity 

The discussion of previous subsection was based on 
qualitative features of (88) which are controlled by the 
IR CFT. To obtain further information regarding proper- 
ties of the Fermi surface and its low energy excitations we 
need to work out (numerically) the data which are con- 
trolled by UV physics. These include: k F (which then 
determines v^ F and 7fc F ), v F , hi and h-2- We will map 
out how they depend on the charge q and dimension A (or 
equivalently mass m) of an operator, i.e. the "phase di- 
agram" 23 of holographic non- Fermi liquids. 

In this section we will often use the notations and equa- 
tions developed in Appendix A. Readers are strongly 
encouraged to read that part first. 

We first solve the Dirac equation (A14) numerically 
with lu = to find k F , for which a/? defined in equa- 
tion is (A38) is zero. This is equivalent to finding the 



bound states of the Dirac equation. See Appendix B for 
more details. For m € [0, i) we consider at the same 
time the alternative quantization of the bulk spinor field, 
whose boundary Green function Gi^m, k) are given by 
Gi.2(— m, — k) (from (A26)). We will thus use negative 
mass to refer to the alternative quantization. The results 
are presented in fig. 5, where we plotted k F dependence 
on charge q for three different masses m = —0.4, 0, 0.4. 
Note in this and all subsequent plots if not stated ex- 
plicitly, we use without loss of generality R = 1, r* = 1, 
g F = 1 and d = 3, for which the chemical potential is 
[i = V3. 

We see in fig. 5 that for a given mass when increasing 
q new branches of Fermi surfaces appear as was observed 
before in [3]. This can be understood from the point 
of view of the Dirac equation as that increasing q al- 
lows more bound states. The lowest bound state has the 
largest kp, 24 which corresponds to the lowest curves in 
fig. 5. We will refer it as the 'primary Fermi surface'. 
For q, m and k large we can also use the WKB approx- 
imation to solve the Dirac equation to find the bound 
states, which is discussed in detail in Appendix B. In 
particular, fig. 12 gives the parameter region where there 
exist Fermi surfaces for large k F ,m,q for various bound- 
ary theory space-time dimensions. 

In fig. 6 we map the values of i>k F , which can be com- 
puted from (56), for the primary Fermi surface in the 
q-m plane. As discussed further in Appendix C, using 
the wave function for the bound state at k F we can eval- 
uate v F and hi,h 2 . Their values for various masses as a 
function of Vk F are presented in fig. 7-fig. 8. 

We now summarize some important properties which 
can be read from these plots and the WKB analysis in 
Appendix B: 

1. For any m > and q, in the standard quantization, 
the existence of Fermi surfaces is always correlated 
with the existence of the oscillatory region which 
requires that for any dimension d (c.f. (67)) 



A < 



\Q\9f 
V2 



d 
2 



(108) 



This is clear from both fig. 5 and fig. 6 for the 
masses plotted there. In fact, the allowed region 
for Fermi surfaces is more stringent than (108). 
This can be seen from fig. 12 (which follows from 
the WKB analysis in Appendix B), which indicates 
that Fermi surfaces only exist for (for any d) 



m 2 R 2 K 



a < 



\q\9F 
V3 



d 
2 



(109) 



This is an abuse of language because the parameters we are vary- 
ing arc couplings in the bulk action, and hence change the bound- 
ary system, not its couplings. 



As discussed in Appendix B, the bound state problem for the 
Dirac equation can be approximately thought of as a Schrodinger 
problem with eigenvalue — fc 2 . 
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At m = 0, in units of the effective chemical poten- 
tial n q , the range for allowed kp is 



d-2 



kp 

= < — < 1 

1) M? 



(110) 



where the lower limit is the boundary of the oscilla- 
tory region. The upper limit in (110) also applies to 
other masses, achieved in the limit m finite, q — > oo. 
The lower limits for other masses are smaller than 
that in (110) as can be seen from fig. 13 (but an 
analytic expression is not known). Note that for m 
approaching the allowed limit (109), kp j \i q lies in 
a small region around 



d-2 



kp 



Hq ^3d{d - 1) 



(111) 



2. For mR G (— |, 0), i.e. for alternative quantization, 
there exists a single Fermi surface which does not 
enter the oscillatory region. This is the primary 
Fermi surface with the largest kp. In fact for the 
small window 



\Q\9f_ 
V2 



< \m\R, 



(112) 



i.e. in terms of boundary theory dimension A 



\m\R, 



d~ 1 A d _ \q\g F 
2 2 V2 



(113) 



there exists a Fermi surface without oscillatory re- 
gion. 

3. For a given m, as one reduces q, kp (and i>k F ) de- 
creases and eventually loses its Fermi surface iden- 
tity by entering into the oscillatory region. Sim- 
ilarly, for a given q, as one increases m, kp (and 
v k F ) decreases and eventually enters the oscillatory 
region. 

4. For a Fermi surface, one can define a topological 
number using the Green's function [21] 



n 



iuj) 



(114) 



which measures the winding of its phase. Here 
Gpi should be considered as a matrix in the spinor 
space and C is any closed loop in the (k± , iuj) space 
around the origin. In our case, for m — since 
detG_R = 1 (see discussion around (A28)), i.e. any 
pole is always canceled by a zero, any Fermi surface 
has a zero winding. For generic m ^ 0, then Fermi 
surfaces here generically have n = ±1. 

5. Except for the single primary Fermi surface for the 
alternative quantization, as indicated in the last 



plot of fig. 5, the Fermi surfaces for the standard 
and alternative quantizations are paired with the 
standard quantization have a larger kp. In the limit 
mR — > i , the paired surfaces for two quantizations 
now have the same kp. 



v F 



From fig. 7, for a given mass and v^ F > \, 
decreases with Vk F . In particular as Vk F — * 5, vf 
approaches zero. From the second line of (92) the 
residue also vanishes in this limit given that hi is 
regular there (see fig. 8). Also note that as v^ F — > 
00, Vp — > 1. 



E. mR 



— 5 limit and free fermions 



We will now consider the mR 



i limit in some de- 
tail (this corresponds to the alternative quantization of 
mR = i), as in this limit the dimension of the operator 
approaches that of a free fermion, i.e. in <i-dimcnsion, 
A — * Note at exactly mR = —i, the bulk wave 

function becomes non-normalizable and there does not 
exist a boundary operator corresponding to it. 25 It is 
natural to ask whether in the limit mR — > — i , the behav- 
ior near the Fermi surface approaches that of a Landau 
Fermi liquid, as recently argued in [4]. 

We first note the following features as mR — > — \ : 



In this limit one finds numerically that for the 
primary Fermi surface, the Fermi momentum ap- 
proaches the upper limit of (110), i.e. 



kp = 



(115) 



2. With kp constrained as in (115), by varying q, v^ F 
can take any values greater than h , 

3. From (C15), one finds that various coefficients 
in (88) behave as 



hi, h-2 oc ( mR 



1 



Vp 



(116) 



This is due to the fact that at mR = — |, the bound 
state wave function becomes non-normalizable with 
a logarithmic divergence. In the limit mR — > — | 
both J 1 and J* defined in (Cll) and (C7) are pro- 
portional to — rrn-- 

With hi — ► and vp — > 1, the Green function (88) 
approaches that of a free relativistic fermion, despite the 
fact that the non-analytic part can still have a nontrivial 
exponent v^ F . Note that equation (115) has a simple 



25 In other words at mR 

rise to a boundary operator of dimension A 



i there is only one quantization giving 



d+l 
2 ' 



18 



interpretation: it is simply the Fermi momentum for a 
free relativistic fermion with Fermi energy fj, q \ 

With hi — > 0, the whole Green function vanishes, sug- 
gesting that at the same time the fermion disappears in 
the limit. Note that nowhere along the limit does a Lan- 
dau Fermi liquid emerge. This picture is consistent with 
general expectations: in the mR — > — ^ limit the mode 
becomes a singleton mode (free fermion) living at the 
boundary and decoupling from everything else. There 
are no bulk degrees of freedom associated with it any- 
more. The fact that we do not see a Landau Fermi liquid 
emerging in the limit is also consistent with our current 
understanding of holography; we do not expect a weakly 
interacting boundary theory to have a bulk description 
in terms of low energy gravity. 

Let us also mention that, at any given mass mR close 
to — 2 j h\ and /12 are small but nonzero. Thus except 
at parametrically small frequencies (i.e. very close to the 
Fermi surface) , the linear analytic term in (88) will domi- 
nate over the non-analytic term. As a result the nontriv- 
ial exponent and the fact that the quasi-particle has a 
finite width will not be easily visible. Turning on a tem- 
perature will generate a new width for the quasi-particle 
and could dominate over the zero-temperature width ex- 
cept at very low temperatures. 

In [4] indications were found that there exists a Fermi 
surface which behaves like a Landau Fermi liquid for mR 
close to at a Fermi energy smaller than us = 0. This 
result is surprising, from the following points of view: (i) 
the Fermi energy is different from (in fact smaller than) 
the value of the effective chemical potential; (ii) it implies 
the existence of some kind of fcrmionic hair outside the 
black hole at us 7^ 0; (iii) from the general philosophy 
of holography mentioned earlier it is surprising to see a 
gravity description of a weakly coupled theory. It would 
be nice to have a better understanding of it. (Note that 
the small w-analysis performed in this paper does not 
apply to any possible Fermi surfaces at a Fermi energy 
different from u> = 0.) 

Finally, we mention in passing that there is another 



limit in which free fermions emerge from the WKB anal- 
ysis. In the limit q — > 00, with mR/q fixed, again one 
finds that /12 — * 0. In this case, depending on the value 
mR/q, vp can take a range of values, see fig. 14. 



F. Double trace deformation and a c-theorem for 
Fermi surfaces? 

As discussed in Appendix A 1, for mR S [0, |) we can 
turn on a double trace deformation O^O in the unstable 
CFT (from the alternative quantization) to flow to the 
stable CFT. We note that the IR CFT of both CFTs 
appear to be the same. It is interesting to examine what 
happens to Fermi surfaces of O under this flow. As an 
example, let us consider m = 0.4. By examining the third 
plot of fig. 5 we observe that in flowing to the stable CFT, 
the primary Fermi surface (which has the biggest radius) 
of the unstable one disappears. As a result, for generic 
q, the number of Fermi surfaces in the stable CFT is one 
smaller than that in the unstable one. For example, at 
q = 2, one starts in the unstable theory with two Fermi 
surfaces with radii given by fci ^> k-z- In the stable CFT 
one finds a single Fermi surface at a radius which is 
greater than but comparable to hi, and much smaller 
than k\. It would be desirable to do a systematic study 
of more examples. 

In the Landau Fermi liquid, Luttinger's theorem says 
that the volume of a Fermi surface gives the number den- 
sity of charge particles. Here it is not clear whether an 
analog of the Luttingcr theorem exists. If we assume 
that the volume of a Fermi surface still has something to 
do with the number of degrees of freedom of the quasi- 
particle excitations, then what we observed in the above 
is consistent with the idea that the number of degrees of 
freedom is always smaller in the IR than the UV. It would 
be interesting to understand this better. It is tempting to 
speculate that some analogue of the c-theorem in terms 
of the volume of Fermi surfaces may exist. 



VII. DISCUSSION 

In this paper we studied the low frequency expansion 
of retarded two-point functions of generic charged scalar 
and spinor operators in a CFT^ at finite charge den- 
sity using its gravity dual, following the earlier numerical 
study of [3]. We showed that the spectral functions ex- 
hibit various emergent critical behavior controlled by an 
infrared CFT described by the AdS2 region of the black 
hole geometry. 



Despite its classical nature, the bulk calculation we 
performed turned out to have an intimate knowledge of 
the quantum statistics of the excitations in the boundary 
theory. In particular, the consistency with the boundary 
theory statistics dictates that a charged scalar field in the 
bulk could have various instabilities including superradi- 
ance, which are absent for a spinor field. We regard this 
as a nice indication of the robustness of the AdS/CFT 
correspondence. We also found a potentially new type of 
scalar instability which appears to be distinct from the 
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FIG. 5: The values of kp as a function of q for the Green function G2 are shown by solid lines for m = —0.4, 0, 0.4. In this plot 
and ones below we use units where R = 1, r* = 1, gp = 1 and d = 3. The oscillatory region, where Uk = -j= \J k' 2 + m 2 — ^- is 

imaginary, is shaded. From (A20) in Appendix A3, Gi(fc) = G2(— k), so &f for Gi can be read from these plots by reflection 
through the vertical k — axis. The m = —0.4 plot corresponds to alternative quantization for m — 0.4 following from 
equation (A26). For convenience we have included in each plot the values of fci? for the alternative quantization using the 
dotted lines. Thus the first (m = —0.4) and the third plot (m = —0.4) in fact contain identical information; they are related by 
taking k — > — k and exchanging dotted and solid lines. Also as discussed after (A29), for m = the alternative quantization is 
equivalent to the original one. This is reflected in the middle plot in the fact that the dotted lines and solid lines are completely 
symmetric. All plots are symmetric with respect to q, k — > —q, —k as a result of equation (A24). 



standard tachyon instabilities (including those induced 
by an electric field). 

Our results suggest a nice description for the low en- 
ergy effective theory near a Fermi surface of a non- Fermi 
liquid. We find at each point on the Fermi surface, there 
lives a CFT 2627 , parameterized by the angle on the Fermi 
surface. This is reminiscent of the Fermi liquid picture, 
except that one replaces the free fermion CFT by a non- 
trivial one. We found a direct relation between the di- 
mension of an operator in the IR CFT and the scaling ex- 
ponent of its spectral function; a relevant operator gives 
rise to unstable quasi-particle excitations at the Fermi 
surface with a zero quasi-particle weight, while an irrel- 
evant operator gives rise to stable quasi-particles at the 
Fermi surface with a nonzero quasi-particle weight. We 
expect this description to be rather general, not restricted 
to theories with a gravity dual, since it involves nothing 
more than general concepts of a fixed point. For example, 
the fact that for marginal operator our expression (97) 
coincides with that of the "marginal Fermi liquid" de- 
scription [8] of the optimally doped cuprates may not be 
an accident and may suggest that the "electron" oper- 
ator is marginal at the possible quantum critical point 
describing a optimally doped cuprates. 

Note that while in this paper we restricted to a charged 
black hole in AdS (which corresponds to CFT at a fi- 
nite charged density), our results should apply to any 
extremal solution with an AdS2 region, e.g. that dual to 
a non-conformal theory at a finite density. One can also 



put the boundary theory at a finite chemical potential 
for some components of angular momenta, whose gravity 
dual is then given by an extremal Kerr-AdS black hole, 
which has an AdS2 region. 

It is worth comparing the form of our Green's functions 
with the well-understood example of a non-Fermi liquid, 
namely Luttingcr liquids in 1 + 1 dimensions, This is in 
some sense a realization of the picture described above, 
in that to construct a Luttingcr liquid one replaces the 
free fermion CFT of each Fermi point with a free bo- 
son of some radius other than the free-fermion radius. It 
differs from our case in two important ways: First, we 
don't know how to obtain the behavior found here as a 
deformation of the Landau theory. Second, the Green 
function we find has a non-analyticity at u> = for ar- 
bitrary k (unlike Gi ut tingcr ~ n^z^ )'t however, only at 
k = kp does this non-analyticity represent a peak in the 
spectral density. 

We now discuss two caveats of our study. The first re- 
gards the stability of this extremal black hole geometry. 
While the black hole is by itself thcrmodynamically and 
perturbatively stable, as discussed earlier holographic su- 
perconductor instabilities 28 can occur if the gravity the- 
ory in which it is embedded contains charged scalars of 
sufficiently large charge or sufficiently small mass [15- 
18]. This might not be of a concern since various phys- 
ical systems to which we might try to apply the mecha- 
nism described here including the normal state of high Tc 
cuprates also exhibit a superconducting instability. Nev- 
ertheless, the criteria for a string vacuum which exhibits 
the fermi surfaces described here but not the supercon- 



As mentioned earlier, the precise nature of the IR CFT is not 
yet understood. It might be a conformal quantum mechanics or 
a chiral sector of a (1 + l)-dimensional CFT. 
In fact, the geometry we are studying displays such a CFT even 
for values of k away from the Fermi surface. 



As mentioned in footnote 14, a neutral scalar with a sufficient 
negative mass square can also condense, whose boundary theory 
interpretation is not yet known. 
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FIG. 6: The 'phase diagram': Shown here are contour plots of the exponent Uk F evaluated on the "primary Fermi surface" 
(the one with the largest fop), as a function of m and q, for each of the two components of the spinor operator. The top and 
bottom rows differ only in the function used to shade the region Vk F £ [0, 5]. Top row is shaded according to a_ while the 
bottom row is shaded according to a+. a± were introduced in (104) and discussed in detail around there. The darker region 
corresponds to larger values of the angles. Both angles are zero at the line of Uk F = \ and increase with a decreasing Vk F ■ 
When an angle exceeds ^ the corresponding pole moves into another Riemann sheet, the regions for which are indicated in the 
plots. As anticipated in the discussion below (104), a_ becomes O(l) only for small q, which a+ becomes small only for vu F 
close to We also indicated the region where there exists an oscillatory region for momentum satisfying (68). This region lies 
above the dashed lines. It is clear from the plots that for m > 0, the region which allows a Fermi surface always lies inside the 
region which allows the oscillatory region. It is also interesting to note that in the left plot the dashed line in fact meets with 
the line for v% F = § at m = — | (not shown in figure). This happens for d = 3 only. For general d dimension, the dashed line 
intersects with m = — \ at Vk F = 2 (d-2) ■ 



ducting instability are reminiscent of those required of a 
string vacuum which describes our universe: one doesn't 
want light scalar fields 29 . In the latter context, a large 
machinery [35] has been developed to meet the stated 
goal, and one can imagine that similar techniques would 
be useful here. It would also be very interesting to un- 
derstand how the condensate affects the Fermi surfaces 
studied here. 

As discussed around equation (17), the black hole solu- 
tion has a finite entropy at zero temperature. Since the 
scmiclassical gravity expression for the entropy is valid 
in the large N limit, this "ground state degeneracy" may 



We thank Eva Silverstein for this analogy. 



be a consequence of the N — > 00 limit. Given that the 
solution is not supcrsymmctric, away from N — 00 (i.e. 
beyond the gravity approximation) likely these states are 
energetically closely-spaced, rather than exactly degener- 
ate. It may be useful to compare this situation to that 
of systems with frustration. While in our geometry this 
non- vanishing ground-state entropy comes together with 
the existence of the AdS2 region which gives rise to the 
IR CFT description, they in principle reflect different as- 
pects of a system and may not correlate with each other. 
For example, suppose the CFTi dual to AdS2 can be 
considered as the right moving sector (with Tr = 0) 
of a 1 + 1-dimensional CFT. Then the nonzero ground 
state entropy should come from the left-moving sector of 
this 1 + 1-dimensional CFT with a nonzero left temper- 
ature. One can certainly imagine a situation where the 
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FIG. 7: The Fermi velocity of the primary Fermi surface of 
various components as a function of vu F > \- Dotted lines 
are for G2. Various values of m are indicated. 



left-moving sector is absent, for which case one will then 
have an IR CFT without a zero-temperature entropy. We 
should caution that it might be hard for such a situation 
to arise as the near horizon limit of a classical gravity 
solution. 

Our results can be generalized in a variety of ways. The 
most immediate is finite temperature, which should shed 
further light on the structure of CFTi. It will be inter- 
esting to follow the instability for scalars to the critical 
temperature [18]. It is also instructive to examine the 
two-point functions of charge current and density fluc- 
tuations. They correspond to fluctuations of the bulk 
gauge field Am in the transverse and longitudinal chan- 
nel respectively. Here we expect even at zero temperature 
there exist hydrodynamic modes like diffusion and sound 
modes as this is a non-neutral plasma and dissipates even 
at zero temperature {e.g. it has a finite entropy). We ex- 
pect our master formula (54) still applies after having di- 
agonalized the perturbations. In particular the diffusion 
and sound modes will depend on the small k expansion 
of a+\ i.e. it should have a zero at k = 0. 

The matching behavior between the inner region and 
outer region is very reminiscent of that for the D3-D7 
system at small densities in [36]. There the geometry is 
more complicated and does not directly give a hint what 
would be the IR theory. Nevertheless it appears likely 
the IR theory there is a non-rclativistic CFT at finite 
density. 
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APPENDIX A: SPINOR CALCULATION 

1. Remarks on dictionary 

A bulk Dirac spinor field ip with charge q is mapped to 
a fermionic operator O in CFT of the same charge. O is a 
Dirac spinor for d odd, and a chiral spinor for d even [33]. 
In both cases the dimension of the boundary spinor O is 
half of that of ip. Since ip has 2^ ~ s ] complex components, 
where [x] denotes the integer part of x, the boundary 
retarded Green function Gr for O is a 2 <l ^~^ x 2[ _ 5 _ 1 
matrix. The conformal dimension A of O is given in 
terms of the mass m of ip by 30 

A = |±mi? (Al) 

where R is the AdS curvature radius. In (Al), one should 
use the + sign for mR > |. For mR e [0, h), there are 
two ways to quantize ip by imposing different boundary 
conditions at the boundary, which corresponds to two 
different CFTs. We will call the CFT in which O has 
dimension A = | + mR the "stable" CFT and the one 
with dimension A = f - mR g (^-, |] the "unstable" 
CFT. In the "unstable" CFT, the double trace operator 
O'O produces a relevant deformation under which the 
theory flows to the "stable" CFT [30]. 

The retarded Green function of O at finite charge den- 
sity can be extracted by solving the Dirac equation for ip 
in the charged AdS black hole geometry. 



2. Dirac equation 

We consider a spinor field in the black hole geome- 
try (2) with a quadratic action 

S= [ d d+1 xy/=g~i(ipr M V M ij-m^) (A2) 



Without loss of generality in this paper we will take m > 0. For 
negative m the discussion is exactly parallel, with m replaced 
by \m\. For odd d the Dirac equation for —m is completely 
equivalent to m as one can change the sign of mass by taking 
i/> — > Tif) where T is the d+ 1-dimensional chirality matrix. For d 
even, different signs of m corresponds to different chirality of O. 
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FIG. 8: hi and h2/\c{kp)\ coefficients in (88) for the primary Fermi surface of various components as a function of v% F . Various 
values of m axe indicated. In the /i2 plot for Gi at m = —0.4: vertical dotted lines denote the minimal v value that is reached 
at q = kp — 0. At this point, our numerical computation becomes less reliable (note that hi appears to blow up there). For 
convenience we also plot |c(fci?)| separately in fig. 9 below. 




0.0 0.5 1.0 1.5 2.0 

v 

FIG. 9: Real and imaginary part of c(&f) as a function of 
Vk F for Gi . The plot for Gi is very similar. Note that due to 
the Gamma function prefactor in (55) the real part diverges 
at half integers. 



with uj a bM the spin connection. Our notations are as fol- 
lows. We will use M and a, b to denote abstract bulk 
spacetime and tangent space indices respectively, and 
(i,i> ■■ ■ to denote indices along the boundary directions, 
i.e. M — (r, n). Underlined indices on Gamma matrices 
always refer to tangent space ones. 
Writing 

V' = (-gg rr )~h- tuJt+tk ' x '^ (A4) 
the corresponding Dirac equation for tp can be written as 

[]tL {Y^d r _ m^T r ) * + iK^T^ = 0, (A5) 

y Qrr 

with 

K fi (r) = (-u(r),k i ) (A6) 



where ip = ip^T- and and u(r) is given by 



V M = d M + \w abM T ah i q A M (A3) « = */ — l^ + Ml-^ll- (AT) 
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As in the case of a charged boson, equation (A5) depends 
on q and (i only through the combination 



fi q = [iq 



(A8) 



which is the effective chemical potential for a field of 
charge q. Similarly, uj should be identified with frequency 
measured away from the the effective chemical poten- 
tial (A8). Due to rotational symmetry in the spatial di- 
rections, we do not lose generality by setting 



ki = k, 



0, ijtl 



(A9) 



Notice that equation (A5) then only depends on three 
Gamma matrices T-, T-, T-. As a result 31 projectors 



n„ 



f 



i - (-l^rq^r 1 ' 



= 1,2, n! + n 2 = i 

(All) 

commute with the Dirac operator of (A5) and 
$„ = n Q *, a = 1,2 



$ = 



(A12) 
i write 

and choose the following basis of Gamma 



decouple from each other. It is then convenient to write 
matrices 



-<7 3 1 

-<7 3 1 

-a 2 l 

a 2 l 



T t - = 



ia x \ 
ia 1 ! 



(A13) 




under which the Dirac equation (A5) becomes 



{d r + 



(ia 2 u+ (-l^fccr 1 ) $ Q . 

/A14) 

In equation (A13), 1 is an identity matrix of size 2~%~ for 
d odd (or size for d even); since the resulting Green's 
functions will also be proportional to such an identity 
matrix, we will suppress them below. Note that (A13) is 
chosen so that equation (A14) is real for real w, k. 

Near the boundary, (A14) has two linearly independent 
solutions given by (with 1/r subleading terms for each 
solution suppressed) 



b a r 



— mR 



a = 1,2 



1J ' 01 \0; 

. (A15) 

To compute the retarded functions, one should impose 
the in-falling boundary condition for $ at the horizon. 



For general ki the projector can be written as 
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i±r^r-fcp-) , 



(A10) 



Then the boundary spinor Green functions have two sets 
of eigenvalues given by 32 



G a {uj,k) = — , a = 1,2 



(A17) 



To see (A17), consider spinors <f>± = |(1 ± with 
definite eigenvalues of T-. Then 



with $ Q 



_ Vol 



2/1 
(A18) 

where we have suppressed zero entries in <p±. Equa- 
tion (A17) then follows using the prescription of [33] (see 
e.g. sec. IIIB). 33 

The dictionary (A17) is for the conventional quantiza- 
tion which applies to any m > 0. For mR G [0, ^), there 
is also an alternative quantization, as discussed at the 
beginning of this section. Similar argument then leads to 



- _ Oa 1 

O a l_r Q 



(A19) 



3. Some properties of the spinor correlators 

In (A14) the equation for $2 is related to that of $1 
by k — > —k, so we immediately conclude that 



(A20) 



As a result the trace and determinant of Gr are invariant 
under k — > — k as should be the case. Given (A20), from 
now on we will focus solely on G±(u),k). For notational 
simplicity, we will also drop the subscript a = 1 below. 
Unless written explicitly all relevant quantities should be 
interpreted as having a subscript a = 1. 

More properties of G can be derived from (A14). For 
this purpose it is convenient to write (A17) as 34 

G = lim e -2mH £L = i, with £=-. (A21) 

From (A14), as in [32] one can then derive a flow equation 
for & 



-2mVM + (it - k) + (u + fc)^ 2 (A22) 



32 That is, when diagonalized, the boundary retarded functions 
have the form 



G R {w,k): 



Gi(w,fc)l 



G 2 0,fc)l 
33 A note on notation: compared to the notation of [3], 



(Ai6) 



where k\ is the unit vector k = k/\k\ 



In (A21) one should extract the finite terms in the limit. 
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with in-falling boundary condition at the horizon given 
by (for uj ^ 0) 



£|r— r 



(A23) 



Properties of G can now be read from those of (A22). 
By taking q — > —q,u> — ► —uj,k — > — fc and £ — > — £ we 
find that the equation for £ goes back to itself, implying 



G(uo,k;q) = -G*(-LJ,-k;-q) 



(A24) 



where the complex conjugation is due to that with u> — > 
—uj the infalling horizon boundary condition turns into 
the outgoing one, which can then be changed back by a 
complex conjugation. 

By dividing both sides of equation (A22) by £ 2 , we 
obtain an identical equation for — | if we also take to — > 



— m, k — > — k. This implies that 
G(w, k; —to) 



G(cj, — fc; m) 



(A25) 



Given equation (A19) we conclude that for alternative 
quantization G can be written as 



G(u), k; m) = G(w, — k; —m) 



(A26) 



That is, alternative quantization can be included by ex- 
tending the mass range for G(uj, k; to) from to > to 



mR > 



Below and in the main text when we speak 



of negative mass it should be understood that it refers to 
the alternative quantization. 

For to = 0, from (A25) and (A20), we find that 



G 2 {u,k) = — 



1 



Gi(uj,ky 

which implies that 

dct G r (uj, k) = 1, 



to = 



TO = 



(A27) 



(A28) 



Note that since a basis change and a Lorentz rotation 
do not change the determinant of Gr, equation (A28) 
applies to any basis of Gamma matrices and any mo- 
mentum. Combining (A20) and (A27) we also conclude 
that at k = 0, 



We will again divide the r-axis into two regions as (37)- 
(38) and consider the low frequency limit (39). The story 
is very much parallel, so we will be brief. 

In the inner region, to leading order in w-expansion 
equation (A30) reduces to equation (D16) of a spinor 
field in AdS2 with to = All the discussion in sec. D 2 
can now be carried over with the replacement v — > v k 



iRl - e 2 ,q 2 - ie, ml = to 2 + 



k 2 R 2 



(A31) 



For example, near the boundary of the inner region, i.e. 

0, the leading order inner solution can be ex- 



t— r* 

panded as 



$? ) (oj : k-a) = v + l * 2 



R 2 ^ Vk 



r — r» 
(A32) 

where v± and can be obtained respectively from (D21) 
and (D22). More explicitly, 



mR 2 ±v k , 
y^Ra+qeJ ' (A33) 



and 



mVk T{-2v k )T{l + v k -iqe d ) 
r(2i/ fc )r(l - v k - iqe d ) 
to + ^ ) R 2 - iqe d - v k 

u 2vk .(A34) 



ikR 



R 2 - iqed + Vk 



As discussed around (D21), Q k depends on the normal- 
izations of v± in (A33). But it can be checked explicitly 
that the final correlation function (54) is independent of 
the normalizations. 

In the outer region we can then choose the two linearly 
independent solutions for the zero-th order equation (i.e. 
(A14) with uj = 0) by the boundary conditions 



(o) 

V± = v ± 



r — r* 



. (A35) 



G 1 (uj,k = 0) = G 2 (uj,k = 0)=i, m = 



(A29) 



Also note that equation (A27) implies that for to = the 
alternative quantization is equivalent to original one [33] . 



The matching, the generalization to higher orders in uj, 
and the low frequency expansion of Gr now work com- 
pletely in parallel as those for scalar fields. More explic- 
itly, perturbatively in to the full outer solution $o can 
be written as 



4. Small-frequency expansion 



$o = r/+ + Gk(uj)r). 



(A36) 



In this subsection, we present the low-frequency ex- 
pansion analysis of section III adapted to the case of a 
spinor field. Now the equation is given by (A14) with 
a = 1 which we copy here for convenience 



(d r + m^g rr (T 3 ^ <J> = . /— — (ia 2 u 

V 9ii 



ka 1 ) <P . (A30) 



with 



r]±\n > 1 arc obtained from solving (A30) perturba- 
tively in the outer region and are uniquely specified by 
requiring that when expanded near r = r*, they do not 



V± 



(o) 
V± 



(i) 



2 (2) 

uj rj y ± > 



(A37) 
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contain any terms proportional to the zeroth order so- 
lutions (A35). Now expanding various r)± ,n > near 
r — > oo as in (A15) 

rf£> r ~°° a^r mR (fj + b^r'™* Q (A38) 

then the retarded function G is again given by the master 
formula (54). 

APPENDIX B: BOUND STATES AT w = AND 
WKB 

1. Scalar bound states 

Setting u> = in (33) defines the outer region differ- 
ential equation for the scalar field. We are then inter- 
ested in examining normalizablc solutions to this equa- 
tion, since this tells us about the spectrum of excitations 
of the boundary theory. This is equivalent to studying 
solutions to this differential equation with boundary con- 
ditions (47) for r]^ and = in (53). The alterna- 



tive quantization window can be achieved by studying 
6+ = 0. As we will sec, bound states will exist for 
a discrete set of momenta, so (at least for the fermion 
problem) these excitations will define a Fermi surface. 

We will study the spinor problem in the next subsec- 
tion. For now the scalar problem will suffice since this 
problem is probably be more intuitive. 

By scaling the wave function and redefining the radial 
coordinate, 4>{r) = ZtJj(s), the wave equation (33) can be 
put in the form, 

-d 2 s i> + V(s)*jj = (-fc 2 )^ (Bl) 

where the "tortoise" coordinate and the rescaling are, 




In doing this we have defined a unique Schrodinger po- 
tential for this problem: 

V = (-u 2 + m 2 g u ) + ((d s In Zf + d 2 s In Z) (B3) 




FIG. 10: The potential for the scalar field defined by (B3) for m 2 = —3/2, —2, —9/4. In terms of the tortoise coordinate s the 
horizon is located at s — — oo and the boundary at s — 0. The oscillatory region is associated with the continuum for s —> — oo 
and "Fermi" surfaces are bound states in the potential well to the right of this continuum. Note the behavior of the potential 
close to the boundary is V(s) ~ (2 + m 2 )/s 2 . 




FIG. 11: A zero-energy Schrodinger potential which is equivalent to the spinor bound state problem (B4). The three plots 
are for m = —1/4,0, +1/4 respectively, for varying values of the parameter k/n q . For m = —1/4 we have rescaled sf — > 
sf/|s»|, Vf — > Vfs 2 so that for all values of fi q /k the potential can be drawn on the same interval —1 < sf < 0. 
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The problem is to find bound states in this potential 
with negative "energy" E = —k 2 . Pictures of this poten- 
tial are shown in Fig. 10 for fixed m 2 and various values 
of q. Examining the r — » r* limit, which corresponds to 
s — > — oo we find that the potential goes to a constant 

V{a) - -^(-(gea) 2 + m 2 R 2 + 1/4) 

Hence there exist a continuum near the horizon for mo- 
mentum satisfying — k 2 > ^ 2 ' R < 2 (— (qed) 2 + m 2 R 2 + 1/4). 
This is the criterion for the oscillatory region. That is, 
we have identified the oscillatory behavior as arising from 
the existence of a continuum. 

Bound states then exist if the potential well in Fig. 10 
close to s — is deep enough. The "sprouting" of bound 
states out of the oscillatory region in the q — k plot of 
Fig. 5 has a nice interpretation in terms of developing 
new bound states as the potential well varies. 



should imagine scaling the potential by dialing k to find 
when a bound states energy eigenvalue crosses zero, this 
will define the fcrmi momentum fcp. In particular for a 
given value of n q /k there are possibly an infinite set of 
bound states energies which cross zero as k is increased. 

Examining the behavior close to the horizon the tor- 
toise coordinate behaves as sp — > — oo for m > and 
sf ^ s, for m < 0. The potential is of the form c/s F 
for m > and c/(sp — s*) 2 for m < where 

(k 2 R 2 R 2 /r 2 + (qe d ) 2 ) , 
C =~ lL 2 R 2 

so that the condition for being in the oscillatory region 
is the usual condition for a singular c/s 2 potential in 
quantum mechanics: c < —1/4. 

3. WKB analysis 



2. Spinor bound states 

The bound state problem is now the first order Dirac 
equation: 



-(m/R)r —k + u 
—k — u (m/R)r 



(B4) 



subject to normalizability conditions for rj^P in (A35) 

and 0$ = in (A38). 

It is harder to have an intuitive grasp over this equa- 
tion, as in the second order problem. We can of course 
square this operator to obtain a Schrodingcr problem, but 
there seems to be no way to define a unique potential, 
such as the one for the scalar. We proceed with an arbi- 
trary choice, in order to give a qualitative understanding 
of the existence of bound states. 

Taking <f> = (y, z) T we can write (B4) as 



\fgllf 

Vgng rr Qd r (Q 



1 d r (Qy) = (u- k)z 



-{k + u)y 



(B5) 
(B6) 



where Q = exp (m J dr '\fg u g rr mr / 'R\ 

write a second order differential equation for tp = z/Q. 
Defining a new tortoise coordinate: 



Then we can 



We can analyze both the scalar field and spinor under 
various limits using WKB analysis. We will focus here 
on the limit q, k, m — > oo with ratios k/q and k/m fixed 
and consider the scalar and spinor problems in parallel. 




FIG. 12: The parameter region where there exists Fermi 
surfaces for large kF,m,q in the WKB approximation, for 
different field theory space-time dimensions d — 3, 4, 6, 50. 
For a given dimension the upper and lower boundary lines 
are defined by the fermi surface moving into the oscilla- 
tory region (lower) and the non-existence of classical or- 
bits (upper.) See Fig. 13 for an alternative way to state 
this in terms of the radii of the two turning points. The 
extreme point to the right of the allowed region is Pd = 

1/V3, (d— 2)/\/3d(d — 1)] and the point at the lower left 



(B7) corner is Q d = (0, (d-2)/ v / d(d- 1) 



we find a zero- energy Schrodingcr equation, 

1 - u/k 



-di F xi> + vh> = 0, (v F )/k z 



1 + u/k 



Q 4 (B8) 



Pictures of this potential are shown in Fig. 11 for 3 dif- 
ferent fixed values of m and various values of n q /k. Note 
in Fig. life has been scaled out of the potential, so one 



As we will see, both problems are governed by the 
"WKB momentum" 



p 2 = k 2 + (m 2 /R 2 )j 



(B10) 



For the scalar field, this is simply the usual potential- 
minus-energy term V — (— k 2 ) from (B3) where terms 
which depend on Z are small so should be dropped. For 
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the spinor it is the negative of the determinant of the 
matrix on the rhs of (B4) . The sign of p 2 will tell us if 
we are in the classically allowed (p 2 < 0) or disallowed 
(p 2 > 0) region. Note that p 2 ~ 0(k 2 ) is large in the 
WKB approximation. In the parameter space shown in 
Fig. 12, there are two turning points ri,r 2 , with a clas- 
sically allowed region in between the two. The WKB 
approximation to the wave function in the three regions 
is: 

• r < ri and r > r 2 (we use a compact notation to 
write both these regions together, with 1,2 corre- 
lated with ± in a self explanatory way) 



V>(r) = 
$(r) = 



^cxpl±| drVff"Srri>(r')J (Bll) 



^1,2 fmr/R^fp 



y/p(k + u) V k + u 
x exp j ± j dr' (^(fg^pir 1 ) + X (r')) 



(B12) 



where \ is an function given by, 

mr'/R 



k + u(r') 
X(r ) = - , dr 



2p(r') 



k + u(r') 



(B13) 



• r\ < r < r 2 



VP 1 



(B14) 



$(r) =- 



D, 



=Re 



mr/R — i 



Vp{k + u) 
where p 2 = —p 2 and 

0B(ri,r) = 
6 F (n,r) = 6 B {ri,r) 



B15) 



dr'Vg u 9rrP(r') 



Mr') 



k + u(r') 



(B16) 
(B17) 



We can formulate a quantization condition by match- 
ing the integration constants (in particular £ = — 7r/4) 
across the turning points using Airy functions. The quan- 
tization conditions for the scalar and spinor problem turn 
out to be, 



7r(n + l/2) = e B {ri,r 2 ) 
7r(n + l/2) = 9 F (r u r 2 ) 



(B18) 
(B19) 



respectively. Note that the quantization condition for the 
spinor only works for m > 0, and no information on the 
alternative quantization region for either spinor or scalar 
can be found with this analysis. 



Using (B19) in Fig.. 13 we plot contours of fixed v 
and q in the WKB parameter region for n = 0. We 
note that the validity of the WKB approximation is for 
n large, however it seems to work remarkably well for 
n = the ground state. One might also imagine that it 
might be exact for n = when the quantization condition 
forces us into the limit kp^q,m — > oo. Such a situation 
does occur at to the upper boundary of Fig. 13 when 
n — ► f2- Actually more care is required in this limit: 
for the scalar field we can formulate a scaling limit in 
which the potential becomes that of a simple harmonic 
oscillator (SHO) located at the radius where r\ — > r 2 . 
In this case, the WKB quantization condition should be 
exact since it is exact for the SHO. Unfortunately it seems 
hard to find the equivalent scaling limit for the spinor, 
see Scc.C3 for more discussion of this. 




FIG. 13: The WKB allowed region for d = 3 showing con- 
tours of fixed v and m (above) and fixed q (below) based on 
the WKB quantization condition (B19) for n = 0. Upper- 
plot: The lowest horizontal contour is v = 0, increasing to- 
wards v — > oo at the upper boundary. The vertical contours 
are for fixed m with m — lying on the kp/fi q axis. The 
contours in this plot demonstrate that the point P3 controls 
the asymptotic slope of the fixed v contours in Fig. 6. Lower 
plot: The contours of fixed q move towards the upper bound- 
ary with increasing q. Note that these contours end on the 
lower boundary (the oscillatory region) . In particular if we fix 
q and increase m, following along the contours in this plot, we 
see that &f decreases until eventually the bound state enters 
the oscillatory region. 
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APPENDIX C: FORMULAS FOR v F ,hi,h 2 
1. Spinors 

In this section we will derive general formulas for vari- 
ous coefficients appearing in equation (88). These coeffi- 
cients are of great importance in characterizing physical 
properties of the Fermi surface as they determine the 
Fermi velocity, the locations of quasi-particle poles, and 
the residues at the poles. The following discussion will be 
very similar to the derivation of the Fcynman-Hcllmann 
theorem, which (not coincidcntally) is commonly used to 
determine dispersion relations in e.g. photonic crys- 
tals [37]. We will focus on the spinor case, and comment 
on the corresponding result for the charged scalar at the 
end. 

Consider the Dirac equation from (A2) 



{T M V M - m)tp = 



(CI) 



in Fourier space where the Dirac operator T> depends on 
lu and fc, which we will collectively denote as A. Now 
suppose that (CI) has a solution ipo for A = Aq, i.e. 



{T M V M \x -m)iPo = 



(C2) 



Consider varying Ao — » Ao + 5X with the corresponding 
solution to (CI) given by i/jq + SXipi. ipi then satisfies 



(r M p M U -m)Vi + r 



M 



&D M 



OX 



^o = 



(C3) 



Multiplying (C3) on the left by dr y/^gipo , integrat- 
ing by parts, and using (C2), we find that 



W(oo) - W(r*) + 
where 



dr 



-gjoT M ^*Po = (C4) 



W 



(C5) 



We will now be more specific, taking A = Ao corre- 
sponding to iq = 0,k = kp and ipo to be that corre- 
sponding to rjjL defined in (A35), i.e. from (A4) ipo = 
(— gg rr )~^r]+ > . Recall that at k = kp, J7+ is normaliz- 

able with a { ° } = 0. We will now consider small oj and k 
variations separately: 

I. Take A = to in (C4). Then we have -01 = 
(—gg rr )~ z T)^L , where r)+ was introduced in (A37), 
and W = ?y+ r-?y+ . Equation (C4) then becomes 



-(0)„r (1) 



iJ% = 



(C6) 



dri-grrg^H^W? ■ (C7) 



where in obtaining (C6) we have used (A38) 

and (A13). Note that near r — ► r*, r)+ ~ {r—r Sf ) Vk , 
thus 



J cx 



1 



2i/ fc -l 



(C8) 



and becomes divergent from integration near r* 
when v>k < \- Since the first term in (C6) is a 
finite constant, the divergence has to be canceled 
by the second term in (C6). As can be checked 
explicitly the second term in (C6) is indeed diver- 
gent 35 for z/fc < | and precisely cancels that of J*. 
For Vk < o the boundary contribution from the 
horizon vanishes. Thus we find that 



(C9) 



where for V). < |, J* should be regularized as dis- 
cussed above. 

2. Take A = k±_ in (C4). Then we have Vi = 

{-gg rr )-id k r]^(k F ), and W = fj^r^d^ . In 
this case the horizon contribution in (C4) vanishes 
and we find 



with 



.7 1 



ibfd k af+ij\, a =Q 



dr V^gipo^ipo 



(CIO) 



(Cll) 



It can be checked explicitly that J 1 is always well- 
defined and finite. We then conclude that 



a (°) J 1 



(C12) 



3. Denoting ip = (—gg rr ) and multiplying 

/ r °° dr^—gipQ on equation (C2) we find that 



(\/^# rr ^o)| r< _ 
which leads to 

a <°> = y~ 



-5^0 r r Vo) 



r 



(C13) 



(C14) 



where v± are given by (A33). 



with 



<;<o 



dry/^gipoT i/j 



Near r - 



(!) /- \ 
rj + ~ (r — r*j 



29 



Using (C9), (C12) and (C14) we thus find that the 
various quantities introduced in (89) can be written as 
(all expressions below are evaluated at k = kp, uj = 0) 



J 1 , (bf? 
v F = Tt , h 1 = — —. 



h 2 = . (C15) 



It can be readily checked that the above expressions do 
not depend the normalizations of v± (even though c(kp) 
does, as discussed after (A34)). 

For v k > \, writing r?^ 0) = M , from (C7) and (Cll) 

we can express vp more explicitly as 

/~<W5^(M 2 -M 2 ) 



Vp 



rdr^grr{-g tt ){\y? + \z\^) 



(C16) 



Since -^—^ = f(r) < 1, the integrands of the numerator 
and denominator pointwise have a ratio less than one, 
from which it follows that vp < 1. This is borne out 
by the numerical results displayed in figure 7. Note that 
the diverging factor of 2v l-\ m equation (C8) causes the 
Fermi velocity to vanish as V). — > i. 



2. Scalars 

Completely parallel analysis can be applied to a scalar. 
One finds that 



•J = q 

r = k F 



drV=g-g u (r,f) 2 



and 



(A - d/2) (&. 



(0)x2 



hi = 

h 2 = \c(kp)\ 



RJ J 1 



(C17) 
(C18) 

(C19) 
(C20) 



and for Vk> \ 



k F Qrygg"^ 

<7 $™drJg{-g*)A t { n f) 



(C21) 



Note that the above equations make it manifest that (for 
q>0) 



h 1 ,h 2 , v F > 



(C22) 



But the value of vp is not obviously bounded. The pos- 
sibility that this 'velocity' may exceed the speed of light 
is not problematic, since the scalar pole represents an 
instability rather than a propagating mode. 



3. Vp,hi,h2 in the WKB approximation 

Here we work in the large q, k,m limit. In particular 
we will take n fixed (and large) such that in this limit the 
two turning points come together at some radius (n — > 
^2) = ^k- We sec this by examining the the quantization 
condition (B19) assuming the two turning points come 
close together. In this case we may approximate p 2 w 
l 2 (r — rx)(r2 — r) with / ~ kp from which (B19) becomes, 



(r 2 - ri) 2 ~ n/k f 



(C23) 



so indeed for fixed n and large kp, the two radii come 
together. Note in Fig. 13 this limit corresponds to the 
upper boundary where the radius rk is determined by the 
ratio mR/qgp. The radius moves from the boundary 
to the horizon as one varies mR/qgp from to 

The wave function will then be localized at this ra- 
dius and expressions for vp,hi, h 2 can easily be derived. 
We will use (C15) to compute these quantities. In the 
expectation values J 1 and J* we may drop the contribu- 
tion from the disallowed regions since it is exponentially 
small. In the allowed region we may replace all oscillating 
functions with their averages, since for large n they are 
highly oscillatory: sin 2 (6 l F - tt/4), cos 2 (#f - 7r/4) — > 1/2 
etc. In particular, using (B15), we make the replacement: 

\ yl 2 ±lz ^D%(u-k F ±(u + kp)) 

p 

Finally the integrals we are left with are over a small 
interval so we may evaluate all smooth functions (such 
as u, g rr , ga . . .) at r = TV The result is: 



J 1 



Di 



J* = —Dp 



2 V grrg ii k F 

2^r 



g Tr g u u 



T2 dr 

— (C25) 

Jr x P 

T2 dr 

— (C26) 

1 P 



The integral in the above equations evaluates to n/l. 
Taking the ratio we find an expression for the Fermi ve- 
locity: 



Vp 



qA t {r k ) 



c rLht( r fc) 



(C27) 



where cii g ht(ffc) = \J f(rk) is the local speed of light at 
radius r^. If we interpret A t {rk) as the local chemical 
potential then this formula is consistent with that of a 
free fermion with a relativistic dispersion relation (with 
the speed of light replaced by the local speed of light). 
Plots of vf are shown in Fig. 14. 

To find h\ and h 2 we need to match the various nor- 
malizations in (B12,B15) to the specified normalization 
at the horizon (A35) with v+ given in (A33). We will 
only be interested in their exponential behavior: 



Dp 
Dp 



C 2 F ~b 



(o) exp(r 2 ) 



c[ ~ cx P (r 1 ) 



(C28) 
(C29) 
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FIG. 14: The Fermi velocity in the limit q,k,F,m large for 
d — 3. The limiting velocity is determined by the radius 
Th (which is in turn determined by mR/qgp) using formula 
(C27). For mR/qgF = the wave function goes towards the 
boundary (r/t — ► oo) and Vf — * 1 ■ This is consistent with the 
asymptotic behavior of vf in Fig. 7, since these plots are for 
fixed m and increasing q. 



so that 



hi ~ cxp(-2r 2 ), h 2 ~ (fc F ) exp(-2r! - 2T C ) 



where the tunneling rates are given by, 



(C30) 



r 2 



lim 



dry/ g u 9rrP - mi?log(r A ) (C31) 



lim / dry/ g u g rr p + v log(r e ) (C32) 



and for h 2 we have included the exponential behavior 



of \c(k F )\ 



-2F C 



Together the factors in (C32) 



should be interpreted as a tunneling amplitude from the 
bound state to the inner region. In the limit considered 
here v is large so eventually the kp 2v term will dominate 
in h 2 which will always be asymptotically small. On the 
other hand T 2 in (C31), once regulated, actually turns 
out to be always negative. Hence h\ will be asymptoti- 
cally large in this limit. 

Unfortunately as it stands it is not clear these results 
apply to the ground state n = with asymptotically 
large k, q, m. However for the scalar field a scaling limit 
exists (similar to the limits considered in the next section) 
where the potential as r\ — > r 2 can be approximated 
by that of a simple harmonic oscillator. By matching 
the SHO wave function onto the WKB wave function 
in the disallowed region, we can show that the answers 
(C27,C30) also apply to any n for the boson. A similar 
result may be derivable for the spinor, but we have not 
been able to formulate it yet. 



APPENDIX D: TWO-POINT FUNCTIONS FOR 
CHARGED FIELDS IN ADS 2 

Here we discuss retarded Green functions for operators 
in a CFT dual to a charged scalar and a spinor field in 
AdS 2 . We will give main results here, leaving derivations 
and more extended discussion elsewhere [191. 



1. Scalars 

We consider the following quadratic scalar action 

S = - J d 2 x^f~g [{D a (j))*D a (j) + m 2 0*0] (Dl) 

with D a — d a — iqA a and the background metric and 
gauge field given by 



ds 2 = ^Hr 2 + dC 2 ) 



(D2) 



Writing 4>(t,() = e iwr 0(w, £), the wave equation for <fr 
can be written as 



-9 C 2 + V(Q4> = 



with 



v i0 = ^-Uf 



(D3) 



(D4) 



To find the conformal dimension of the operator O dual 
to 0, we solve (D3) near the boundary ( — > and find 
that 

= A(i-» (1 + 0(C)) + BC^ (1 + O(0) , C -> 

(D5) 

with 



m 2 R\ — q 2 e 2 d 



1 



(D6) 



Since u> in (D3)-(D4) can be scaled away from re- 
defining C, we conclude that in (D5), A ~ and 
B r~j u)2+ u and thus (after imposing infalling boundary 
condition on <p at the horizon) 



g R (u) oc — uj 2v 



(D7) 



which implies a coordinate space correlation function 



with the conformal dimension 5 of O given by 
. 1 



:ifi 



(D8) 



(D9) 



36 One can also reach the same conclusion by assuming a boundary 
coupling J dr <j>oO and considering a conformal scaling in the 
boundary theory. 
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Notice that dimension 5 also depends on charge q. In 
particular, it is possible for v to become imaginary when 
q is sufficiently large. Physically this reflects the fact 
that in the constant electric field (D2) particles with a 
sufficiently large charge q can be pair produced. 3 When 
v is imaginary, there is an ambiguity in specifying Qr 
since one can in principle choose either term in (D5) as 
the source term. We will follow the prescription as deter- 
mined by the —ie term in (D6), as this will be the choice 
one needs to use when patching the AdS2 region to the 
outer region of the full black hole geometry. 

Equation (D3) can in fact be solved exactly and one 
finds the full retarded Green function is given by [19] 



Gr{uj) = 2ve~ 



Y{-2v)Y(\+v-iqe d ) 



(2w) 



2,i 



Y(2v)Y (\-v- iqe d ) 

(D10) 

Equation (D10) has the form of the retarded two-point 
function of a scalar operator in a (l+l)-dimensional CFT 
with left/right-moving dimensions and momenta 



Sl = 8r = - + u, 



PL = q, 



PR = w 



(Dll) 



in a (l+l)-dimcnsional CFT with left/right temperatures 
given by 



Aired' 



T R = . 



(D12) 



Thus it is tempting to interpret the CFTi dual to AdS2 
as the right moving sector of a (l + l)-dimcnsional CFT. 38 
Also note that the advanced function is given by 



Y{2v)Y(\-v + iqe d ) 



(D13) 



2. Spinors 

We consider the following quadratic action for a spinor 
field ip in the geometry (D2) 

S = J d 2 x^gi(ipr a D a i}- m^tp + im't/jTip) (D14) 

where we have included a time-reversal violating mass 
term proportional to rh which in our application will be 
related to momentum in R d_1 . 

It is convenient to choose the following Gamma matri- 



ces 



:!') 



ri = 



(D15) 



where underlined indices again denote those in the tangct 
frame and a 1 are standard sigma matrices. The equations 
of motion for ijj can be written in Fourier space as 



ia 2 + (ma 3 + rha 1 ) $ (D16) 



with $ = (— gg^) Near the boundary £ — > 0, equa- 
tion (D16) becomes 



C<9 C $ = -17$, U 



mR 2 mR 2 - qed 
rhR 2 + qed -mR 2 



(D17) 



Thus as £ — * 0, $ can be written as 

$ = Av+C v (1 + 0(C)) + Bv-C (1 + O(0) (D18) 

where v± are real eigenvectors of U with eigenvalues ±v 
respectively and 



v = J (to 2 + m 2 )R. 



q 2 e 2 d 



(D19) 



Imposing the infalling boundary condition for $ at the 
horizon, the retarded Green function for the boundary 
operator in the CFTi dual to ip can then be written as 40 



(D20) 



again suggesting the operator dimension to be given by 



5 = 



v. There is, however, an ambiguity in (D20) 



as the ratio depends on the relative normalization of v± ; 
if we take v± —* \±v±, then Qr — ► j^-Gr. This ambi- 
guity 41 will not be relevant for the present paper as we 
will see later the ambiguity cancels in the matching pro- 
cedure and the correlation function for the full geometry 
will not depend on the normalization. 

As in the scalar case differential equation (D16) can be 
solved exactly and with the choice of v± by 



v± 



TO-/? 2 i V 



iflRj 



qed 



(D21) 



we find that Qr can be written as [19] 

_ im/ Y{-2v)Y{l + u-iqe d ) 



Y_{2v)Y{l-v-iqe d ) 

(m + irh) R? — iqed — v, K o„ 
: - — — „ — (2uj) 2v . (D22) 



(m + irh) R 2 — iqe d + v 
The advanced function is given by 

, v Y{-2v)Y{l + v + iqe d ) 



Ga{u) 



Y{2v) Y (1 - v + iqe d ) 



As we discuss in the main text, when embedded in the full theory 
this causes an instability for scalars, but not for spinors. 
We should caution, however, that many other aspects of this 
theory should be studied before one can really draw a conclusion. 
They are chosen to be compatible with the choice made in ap- 
pendix B, with (D2) arising as the near horizon limit. 



As discussed earlier for a bulk spinor field the number of compo- 
nents of the boundary operator is always half of that of the bulk 
field. 

It will be discussed in more detail in [19]. 
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(m - irh) R 2 + ige d - v (C) 
(m — irh) R 2 + iqed + v 



(2w) 2 " . (D23) 



Equation (D22) is again suggestive of spin-i operator 
with lcft/right-moving dimensions and momenta 



l + i/, 



OR = ~- s rV, PL=q, VR = U 

(D24) 

in a (1 + l)-dimcnsional CFT with temperatures given 
by(D12). 

Finally note as with scalar case, (D19) can become 
imaginary when q is sufficiently large, in which case the 
prescription for determining Gr is again determined by 
the —ie term. 



For spinor from (D22) and (D23) we find 



Gr(w) _ 27TW cos it (is + iqe d ) _ e" 



3 -27rge ti 



Qa[oj) cos tt(v - iqe d ) e 2mv + e" 27 ^ 

(D26) 

Writing Qr = cuj 2 ", then for real is equations (D25) 
and (D26) give the phase of c and for imaginary v they 
give the modulus of c. 



A useful formula 



4. Finite temperature generalization 



Here we give a nice formula for Qr which can be used to 
derive equations (60)-(63) in the main text. For scalar, 
using (D10) and (D13) we find that 



Ga{uj) 



-27tiv 



siTtiriy + iqed) 
sin7r(V — iqed) 



-2irqe d 



^2-rziv 



-2nqe d 

(D25) 
I 



One can in fact generalize the above discussion to finite 
temperature, i.e. to the AdS2 part of the metric (24). 
Details will be given in [19]; we mention the result here 
because it gives useful information about the analytic 
structure of the Green's functions in u at zero tempera- 
ture. One finds respectively for a scalar and a spinor 



Gr(lo) = (2ttT) 



2 *r(-2i/)r 



+ iqe d ) r (i + v - iqe A ) 



r(2^)r (i - v - ^ + iqe d ) T(\-v- iqe d ) 



(D27) 



G R (w) = (2ttT) 



2v r(-2i/) r(± + v - + iqe d ) r (1 + v - iqe d ) (m + irh)R 2 - iqe d 



T{2v) r 



2^ + iqed) T(l-u- iqe d ) (m + ifh)R 2 - iqe d 
I 



(D28) 



Note that the branch point at lu = of zero temperature 
now disappears and the branch cut is replaced at finite 
temperature by a line of poles parallel to the negative 



imaginary axis. In the zero temperature limit the pole 
line becomes a branch cut. Similar phenomena have been 
observed previously [38] . 
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